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IN HONOUR OF THE SCIENTIST 





21 June, 2020, professor Rudolf Anatolievich Neidorf, a remarkable scientist, 
engineer and organizer, Dr.Sci. (Engineering), died unexpectedly at the age of 76. He 
would have turned 76 years old 17 October. 

Rudolf Anatolievich Neidorf was an expert in the field of automatic and 
automated systems, well-known in Russia and abroad. R.A. Neidorf embarked upon his 
scientific career at Novocherkassk Polytechnic Institute (SRSTU, Novocherkassk) as an 
expert in the automation of chemical technology processes, and he continued his activities 
at Rostov State Agricultural Academy and Don State Technical University (DSTU). 


R.A. Neidorf’'s main area of expertise was the automatic control theory, in 





particular, development of the methodology for the structure synthesis of nonlinear 
automatic control systems. This is a complicated research area at the intersection of intuition, art and mathematics 
which requires nonstandard ideas and scientific courage to stand up for them. R.A. Neidorf’s PhD and habilitational 
theses, monographs, scientific publications and scholary papers, as well as the dissertations of many of his disciples, are 
devoted to the research in this area. A wide range of scientific interests of R. A. Neidorf included also tasks and methods 
of mathematical modeling, optimization, and information processing. R.A. Neidorf developed these problems together 
with his graduate students. The joint research results were recorded in the publications and manuals. 

The results of the theoretical and R&D work of Rudolf Anatolievich were tested at the major Russian and 
foreign scientific conferences and introduced into production. One of the examples is the method of experimental 
building of mathematical models of highly nonlinear systems proposed by R. A. Neidorf and a group of young 
scientists. The solution found application in the selection of structural elements of the MAAT transport system 
(Multibody Advanced Airship for Transport) (2012—2015). It was the basis for a scientific project on the theoretical 
underpinning of the method of experimental building of mathematical models and its implementation in scientific 
practice under a grant from the Russian Foundation of Fundamental Research. 

R.A. Neidorf made a significant contribution to the development of engineering education in Russia, the 
improvement of the educational process, the training of engineering and scientific top-qualification brainpower. Two 
Dr.Sci. and more than twenty Cand.Sci. in the field of mathematical modeling, control theory, and operation research 
methods, trained by Rudolf Anatolievich, work at the universities and enterprises of the Rostov region. 

R.A. Neidorf organized the training of engineering personnel majoring in “Management and IT in Engineering 
Systems”, scientific brainpower majoring in “Information Systems Software and Administration” and “Software 
Engineering” in DSTU. 

In 2002—2005 and 2006-2013, R.A. Neidorf was Head of the “Information and Control Systems” Department, 
Rostov State Agricultural Academy, and “Computer and Automated Systems Software” Department, DSTU. He made a 
great contribution to the scientific activities of the higher educational institutions, which later merged into a regional 
basic university, DSTU, and intensified R&D in the field of information control problems. R.A. Neidorf, in particular, 
developed optimization methods, control laws for various engineering systems, as well as search heuristic problem 
algorithms; he established the priority of R&D in this area. 

At that time, he forged close links with the “Engineering Systems Control” Department, Kislovodsk 
Humanitarian Technical Institute, and the “Automated Information Systems” Department, Ukhta State Technical 
University, and organized visiting professorship. Lately, he had done a lot of work at the Center for Scientific 


Competence while developing methods for training top-qualification brainpower. 
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Since 1996, R.A. Neidorf headed the specialized Dissertation Board D.212.058.02. In 2010, he organized the 
annual seminar-conference “System Analysis, Management and Information Processing”. Rudolf Anatolievich was a 
member of two editorial boards: “Vestnik of Don State Technical University” and the collection of research papers of 
the seminar-conference organized by him. R.A. Neidorf was author and co-author of more than 500 scientific 
publications including 5 monographs (3 foreign ones), 12 textbooks and teaching aids, more than 70 patents and 
copyright certificates. 


Not many people know that Rudolf Anatolievich was a poet and editor of the monthly Ataman News, author of 
158 poetical compositions. 
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Transverse vibrations of a circular bimorph with piezoelectric and 
piezomagnetic layers 





A. N. Solov'ev, Do Thanh Binh, O. N. Lesnyak 


Don State Technical University (Rostov-on-Don, Russian Federation) 


Introduction. Transverse axisymmetric oscillations of a bimorph with two piezo-active layers, piezoelectric and 
piezomagnetic, are studied. This element can be applied in an energy storage device which is in an alternating magnetic 
field. The work objective is to study the dependence of resonance and antiresonance frequencies, and electromechanical 
coupling factor, on the geometric parameters of the element. 

Materials and Methods. A mathematical model of the piezoelement action 1s a boundary value problem of linear 
magneto-electro-elasticity. The element consists of three layers: two piezo-active layers (PZT-4 and CoFe,Ox,) and a 
centre dead layer made of steel. The finite element method implemented in the ANSYS package is used as a method for 
solving a boundary value problem. 

Results. A finite element model of a piezoelement in the ANSYS package is developed. Problems of determining the 
natural frequencies of resonance and antiresonance are solved. Graphic dependences of these frequencies and the 
electromechanical coupling factor on the device geometrics, the thickness and radius of the piezo-active layers, are 
constructed. 

Discussion and Conclusions. The results obtained can be used under designing the working element of the energy 
storage device due to the action of an alternating magnetic field. The constructed dependences of the eigenfrequencies 
of the resonance and antiresonance on the geometric parameters of the piezoelement provide selecting the sizes of the 
piezo-active layers for a given working frequency with the highest electromechanical coupling factor. 


Keywords: energy storage device, piezoelectrics, piezomagnetics, finite element method, natural oscillation frequencies. 
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Introduction. In the development of sensor and measuring systems, modern small-sized household appliances, 
cellphones and wireless sensor systems, powerful sources of energy are not required for monitoring and diagnosing the 
technical condition of various objects, but mobility and nonvolatility of the above devices are mandatory. Energy 
storage devices with piezo-active elements that directly convert the energy of mechanical vibrations into electrical 
energy are widely used to power this kind of apparatus. In [1-3], energy storage devices using piezoelectric generators 
under the action of mechanical loads are studied. 

If the system is in an alternating magnetic field created by permanent magnets mounted on rotating parts of the 
machine, then the piezomagnetic layer is deformed along with the piezoelectric element. Due to this, an electric current 
is generated. In [7,8], Guo-Liang Yu and others discussed theoretical models of multilayer magnetoelectric composites 
for the magnetoelectric response at resonant frequencies corresponding to vibrations of the bending and extensional 
modes. A theoretical model of magnetic energy collectors using a functionally gradient composite cantilever was 
analyzed to improve the collecting ability and adjust the resonance frequency in [9]. In [10], Y.F. Zhanga and others 


studied bifurcations, periodic and chaotic dynamics of a tetrahedral composite multilayer piezoelectric rectangular plate 
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with simple supports. In [11], bending and free vibrations of a magnetoelectroelastic plate with surface effects were 
studied. 

In this paper, axisymmetric vibrations of the device are considered. The effect of the radius and thickness of 
the piezoelectric and piezomagnetic plates on the frequency characteristics of the device (eigenfrequencies of resonance 
and antiresonance) and on the efficiency of the conversion of vibratory energy into electrical energy, which is 
characterized by an electromechanical coupling factor, is studied. As a method for solving the problem, the finite 
element method implemented in the ANSYS package is selected. 


Z - axis of Symmetry 
a 

















hp OW Piezoelectric plate 
u,=0 
T=0 


i Piezomagnetic plate 
Metal plate 


Fig. 1. Scheme of half axial section of piezoelectric generator 


The energy storage device under consideration is an axisymmetric design, which consists of a metal disk 
(substrate) with two piezoelectric layers glued on it (Fig. 1). The top layer is piezoelectric; the bottom one is 
piezomagnetic. The flat surfaces of the piezoelectric layer are covered with electrodes that are connected to an external 
electrical circuit, or one of the electrodes is free while the other is set to zero electric potential. The piezomagnetic layer 
is affected by an alternating magnetic field according to the harmonic law; the outer radius of the substrate is pivotally 
fixed. The mathematical model of transverse steady-state vibrations of the described construction is the boundary-value 
problem of the linear theory of piezo-magnetoelectric elasticity [4]. The ANSYS package, which implements a 
piezoelectric model, is selected as a solution tool. 

The boundary value problem for a piezo-magnetoelectric body consists of a system of equations and boundary 
conditions [4]: 


V-o+pf=pu, V-D=0,, V-B=0 (1) 
o=c:e—-e -E-h’-H 
D=e:6e+K-E+a-H (2) 


B=h:e+o' -E+p-H 
e=—(Vu+(Vu)'), E=-Vo, B=-V6o, (3) 
where o and é€ are tensors of mechanical stress and strain; D and E are vectors of electric induction and electric field 
strength; B and H are vectors of magnetic induction and magnetic field strength; p is density of the material; c is 
tensor of elastic modules; e is tensor of piezoelectric modules; h is tensor of piezomagnetic modules; « is tensor of 
dielectric permeabilities; a is tensor of magnetoelectric modules; uw is permeability tensor; f is vector of the density 
of mass forces; 6, is volume density of electric charges; u is displacement vector; ~ and 9 are electric and magnetic 
potentials. 
Boundary conditions are specified for the mechanical, electrical, and magnetic components of the fields. 
Mechanical boundary conditions. Let the surface S consist of two parts T, andI’,, so thatS =T, UTl,, 
moreover, AT, =2. 
u=U onTI,, n-o=p onl,. (4) 
Electrical boundary conditions. Let the surface S consist of two parts T, andI,, so thatS=T,UT,, 
moreoverl, 01, =. 
p= (x,t) on T;, n-D=—o, on I,, (5) 
where o, is surface-charge density. In addition, if the electrodes are connected to an external circuit, two conditions 
should be added: 
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Ol, =v, [[n-DdS=I, (6) 


Sr 
where s, is electrode area; v is unknown potential which is found from the second condition; / is current. 
Magnetic boundary conditions. Let the surface S consist of two parts I; andI,, so thatS =I; UD,, moreover 
rPeor., =O. 
d=(x,t) on T,;, n-B=o, on I, (7) 
where o, is density of free surface currents along the boundary. 
For the elastic layer, the first equations in the system (1)-(3) are used; the components of the displacement 


vector u are unknown. For the piezoelectric layer, the second equations are added to them; u and the electric potential 


© are unknown. For the piezomagnetic layer, the third equations are added to the first equations, and the magnetic 
potential is added to the unknown displacements. In this case, the relations (2) are transformed through zeroing the 


corresponding constants. 

Finite element modeling. The computer model of the device is built in the ANSYS finite element package. 
The metal substrate (steel) has thickness fA and radius r. The piezo-active layers consist of one piezoceramic and one 
piezomagnetic plates which are polarized in thickness, have thickness h, and radius r, (Fig. 1). 

In the finite element model, devices for the metal and piezoelectric layers are used as finite elements PLANE42 
and PLANE13, respectively. In this work, the piezomagnetic layer is modeled by the finite element PLANE13, in which 
the piezoelectric properties of the material are replaced by piezomagnetic ones. This can be done for two reasons: the 
piezoelectric layers do not contact each other; qualitatively, the equations for the electric and magnetic potentials 
coincide. 

The properties of the piezoelectric layers used in the calculations are presented in Tables 1—2: piezoceramics 
— PZT-4, piezomagnetic material — CoFe204 [5, 6], adhesive layers are not taken into account. 

Table 1 


Material properties of PZT-4 piezoceramics 


e e e 
C",Hpa | C*%,Hpa | C..,Hpa | Ci, Hpa an al i | ome: | aw a 


77.8 74.3 iis 25.6 730 635 


€, = 8.85x10 F/m, PZT-4 piezoceramic density p = 7500 kg/m’ 





Table 2 
Material properties of CoFe,O, piezomagnetic element 


Ci. Cae Gas oe Cra Q;,, N/A Q;;, N/A O;;, N/A Mir» M3 > 


H Hpa N s/KI’ Ns/KI 


Hpa Hpa Hp 
2 


a pa m m m 
286 269.5 | 45.3 | 580.3 699.7 5.9x10" a0 





CoFe,O, piezomagnetic element density is p = 5290 kg/m° 
The elastic properties of the isotropic substrate material are characterized by the Young's modulus F and the 


Poisson's ratio v, E = 200 Hpa, v = 0.29 , density p = 7860 kg/m’ (steel) were used in the calculations. 


To achieve high accuracy in the calculations, the size of the final element of the metal layer is set to the value 


not higher than 1/5 its thickness; the size of the final element of the piezoelectric layers is automatically set. The finite 


element grid of the energy storage device 1s shown in Fig. 2. 
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Fig. 2. The finite element grid of the energy storage device. 
Surfaces 1, 4 are outer surfaces of piezoelectric layers; 
surfaces 2, 3 are contact surfaces between piezoelectric layers and metal layer. 


When analyzing the natural vibrations of a piezo-electromagnetic energy storage device, it is assumed that the 
following mechanical and electromagnetic boundary conditions are met. At the left end, symmetry conditions are 
specified, the right end is pivotally fixed (Fig. 1). To calculate the resonance frequency of the device, electric potentials 
are set on surfaces 1, 2; on surface 3, the magnetic potential is set; the magnetic-flux density is set on surface 4. In the 
case of calculating the antiresonance frequency of the device, the same boundary conditions as for calculating the 
resonance frequency are set on surfaces 1-4. However, the electric potential on surface 1 is unknown and is found from 
the condition (6). 

Numerical results. The natural vibrations of a piezoelectric element are considered, the radius of which is 
r,—9.8 mm, the thickness of the piezoelectric layers is h,=0.5 mm, the radius of the substrate is r=10 mm, and the 
thickness of the substrate is h=0.1 mm. Table 3 presents the first three natural frequencies of antiresonance. 


Table 3 
Antiresonance natural frequencies 


Antiresonance 


Distribution of vertical displacement on deformed state and 
distribution of electric and magnetic potential on undeformed sample 


frequency 
kHz 
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Antiresonance 


Distribution of vertical displacement on deformed state and 
distribution of electric and magnetic potential on undeformed sample 


frequency 
kHz 


115,372 





Table 3 presents the distribution of vertical displacement, electric and magnetic potential while the latter are 
presented in a single scale, here, the display of the magnetic potential is not visual. To build an approximate theory of 
calculating the vibrations of the considered bimorph, an analysis of the stress-strain state, electric and magnetic fields, 
which shows that for lower bending modes, there can be accepted hypotheses about their distribution corresponding to 
the bending of the plates. 

In the results of numerical calculations presented below, the dependence of the eigenfrequencies of resonance 
and antiresonance, and the electromechanical coupling factor on geometric parameters is investigated. 

When the number of alternating piezoelectric and piezomagnetic layers becomes sufficiently large, it 1s 
possible to use an approach based on the effective properties of the piezo-magnetoelectric composite [5, 6]. In this case, 
all equations of the problem (1) - (3) are used. 

The thickness value of the piezoelectric layers h, varies in the range of 0.3 + 0.7 mm, and the radius value r, — 


in the range of 6.8 + 9.8 mm. 


fal 





Fig. 3. Natural frequencies: a) resonance frequency, 
b) antiresonance frequency 


Fig. 3 shows the natural frequency dependences on the thickness h, and the radius r, of the piezoelectric 
layers. Fig. 3 shows that the values of the natural frequencies increase with increasing the radius. 

Fig. 4 shows the dependence of the electromechanical coupling factor on the thickness h, and the radius r, of 
the piezoelectric layers. Fig. 4 shows that the value of the natural frequency increases with increasing the radius of the 
piezoelectric layers r,, but decreases with increasing thickness h, of the piezoelectric layers. 
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Fig. 4. Electromechanical coupling factor 


Conclusion. An axisymmetric finite element model of an energy storage device based on round plates in the 
ANSYS package is considered. The active elements of the energy storage device are the piezoelectric and 
piezomagnetic plates fixed on a metal plate. The effect of the geometric characteristics of the piezoelectric layers under 
certain restrictions on the device dimensions, on the electromechanical coupling factor, which shows the efficiency of 
the energy storage device, is studied. The calculation results presented in the paper allow selecting rational sizes of 
piezoelectric elements operating at a certain frequency and having the greatest efficiency. This analysis shows that the 
maximum value of the electromechanical coupling factor is achieved when the thickness and radius of the piezoelectric 
layers take the largest and smallest values, respectively, within the considered limits. 
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Investigation of crack propagation in the surface white layer of rail steel ut 10 


A. Yu. Perelygina, V. Yu. Konyukhovy, A. E. Balanovskii 
Irkutsk National Research Technical University (Irkutsk, Russian Federation) 


Introduction. The paper is devoted to the evaluation of cracking of white layers formed on the surface of the rail while 
in operation. Cracks are detected in the white layer of rail steel after one thousand test cycles. This is due to tensile and 
shear stresses on the surface of the wheel—rail contact spot. The paper presents the study results of the morphological 
characteristics of the white layer on the rail surface. 

Materials and Methods. The object of study (rail surface after operation) was examined under a microscope. Then, a 
two-dimensional model of finite elements of the plane deformation was developed to simulate the dynamic 
characteristics of the white layer cracking. Mathematical models describing crack propagation are proposed. For this, 
we applied the criterion of the elastic plastic fracture mechanics, the J-integral method. The SYSWELD program 
performed numerical modeling of the formation of a white layer and the distribution of residual stresses. 

Results. Optical images of the microstructure of the cross section of a white layer on the rail surface after operation are 
presented. Two different types of cracks were fixed at the trailing edge of the white layer of the samples studied. The 
SYSWELD program visualized fragments of simulating the mechanism of the white layer formation with the distribution 
of residual stresses, compression, and tension. The calculation results show that the values of the J-integral for all three 
cracks slightly decrease if the crack length reaches 10—50 um. 

Discussion and Conclusions. The results obtained are applicable to assess the wear resistance of rail steels and predict 
the direction of crack growth. Comparisons of J-integral maxima have shown that under identical load conditions, crack 
no. | is likely to grow faster than cracks nos. 2 and 3. With an increase in the length of the crack, the maxima of the J- 


integral of all three cracks decreased. 
Keywords: rail steel, white layer, crack, elastic plastic fracture, J-integral method, distribution of residual stresses. 
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Introduction. It is known that parts of machines and mechanisms, various functional structures often 
experience catastrophic brittle fracture. Moreover, depending on the operating conditions, the metal may undergo 
plastic or brittle fracture [1-3]. Brittle fracture occurs due to the growth of cracks that suddenly become unstable and 
propagate in the material at sonic velocity. Cracks in the metal can be of technological origin or appear and grow during 
operation. White non-etching layer (WEL) is a phenomenon that occurs on the surface of an operated rail under the 
impact of wheels. WEL is formed due to strong multicycle plastic deformation. It has been established that due to 
fatigue from contact rolling on the rail surface, cracks bind to the white layer [4-6]. The layer is named because of its 
resistance to acid etching during metallographic preparation and the white “faceless” appearance (it looks like this 
under a microscope). A white layer is usually found in contact spots on the rail surface. Its depth is ~ 0.10—100 um after 
several months of rail operation [4]. An important feature of the white layer is hardness. It reaches 1300 HV [5], but is 
usually in the range of 700-1200 HV [6]. This condition can cause the formation of brittle cracks in the white layer and 


subsequent propagation of fatigue cracks. The authors [4], having examined the rails after operation, found that cracks 
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were present in the white layer, but did not penetrate the interface with the base metal of the rail [4, 5]. Cracks 
developed on the surface and propagated along the interface between the WEL and the pearlite structure of the material 
causing severe wear. The mechanism of the formation of the white layer and its microstructure are studied by many 
authors [6-12]. 

The microstructure of the layer 1s identified as the following mixtures: 

- martensite with residual austenite and pure martensite saturated with nitrogen and carbon; 

- ferrite with martensite; 

- residual austenite and martensite with tempered martensite; 

- ferrite with cementite, carbide and nanocrystalline phase a-Fe. 

It was found [6, 7, 9] that the white layer has a martensitic microstructure with a high dislocation density. The 
authors of [8—10] suggested that WEL consists of a nanocrystalline phase of a-Fe with grain sizes from 15 to 500 nm. 
According to [7], WEL consists of strongly deformed perlite, nanocrystalline martensite, austenite, and cementite. In 
[8], the results of X-ray measurements of residual compressive stresses in WEL are presented, and significant 
compressive stresses of the rail (~ 600 MPa) in both directions, longitudinal and transverse, are specified. An analysis 
of the papers [4-19] shows that cracks in the white layer of rail steel are detected after all tests, even those that were 
completed far just one thousand cycles. The reason for this is mixed tensile and shear stresses on the surface of the 
wheel — rail contact spot. After the appearance of a crack in the white layer, it grows rapidly due to the fragility of the 
martensite structure until it reaches the interface with the base metal of the rail having a perlite structure — here a 
different microstructure deflects the crack and slows its propagation. This deviation is caused by the orientation of the 
cementite plates in the pearlite structure [10—19], which are parallel to the rolling surface. According to [12], a decrease 
in rail fatigue resistance is associated with contact pressure and slip coefficient. The authors of [13] have found that the 
initiation of a crack in a white layer is caused by the traction force and the designed transverse shear load to the rail 
track surface. Two types of cracks are identified: 

- leading defect crack in the contact spot (caused by a shift); 

- back crack (brittle fracture in the form of a wedge). 

The authors of [14] associated the short-wave dynamic interaction between the wheel and the rail with the 
beginning and growth of surface deformation. Under the conditions of contact rolling in rail steel under fatigue, three 
types of cracks were identified [15]: 

1) perpendicular contact surfaces, 

2) with a tilt angle, 

3) parallel contact surfaces at various depths. 

To determine the stress state near the crack tip under conditions of contact with rolling, the authors proposed a 
numerical method and recorded a change in the crack shape. In [16], the criteria for the direction of crack growth in the 
mixed mode were used to evaluate cracking behavior. It is shown that cracks grew in the direction of the plane of 
maximum shear stress, and not perpendicular to the plane. The authors of [17] developed a two-dimensional finite 
element model for simulating the crack behavior in the rolling contact (wheel and rail) with four short cracks. It was 
found that the shear stress plays a dominant role in the crack growth, a longer crack grows upward, and this causes a 
breakaway of the surface layer. Using a two-dimensional computational model, the authors of [18] found that the 
growth rate of cracks on the rail surface increases with increasing crack length and starts to decrease after a certain 
depth. In the study on the dynamic wheel and rail interaction [19], a three-dimensional model of finite elements was 
used and various cracking behavior was recorded due to the static and dynamic solutions. In [20], the residual stress 
effects were taken into account when studying the path and rate of crack growth. Many works including [21, 22] show 
the influence of a liquid on the further development of cracks. It is assumed that the liquid on the surfaces of the cracks 


significantly affects their opening (stress intensity factor of mode I). If the crack develops under dry conditions, then the 
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shear regime dominates [22]. Until now, linear mechanics of elastic fracture has been widely used under the numerical 
simulation of the contact of wheels and rails [23, 24]. This approach evaluates the resistance to unstable crack 
propagation, or the crack resistance of metals, under static loading according to one or more fracture criteria: 

- force (critical stress intensity factor K ,,); 

- deformational (critical opening at the crack tip); 

- energy (the critical value of the J-integral, the work of plastic deformation and fracture). 

In practice, materials near the contact area can be plastically deformed due to the high axial load [3—6]. To obtain 
a more accurate solution, calculations using elastoplastic fracture mechanics or energy-based criteria should be used. In 
addition, cracks can be formed in areas adjacent to the white layer [7—12], and the mechanisms of their development are 
still not properly understood. Papers on the formation of a white layer and cracks [14, 15] are published; however, the 
behavior of the crack propagation in the area of the white layer is not studied in detail. Moreover, in relation to the 
situation under consideration, the impact of loading conditions, friction, and other interaction parameters is not 
described. This paper presents results of the study on the morphological characteristics of the white layer on the rail 
surface. 

Materials and Methods. At first, the object of study was examined under a microscope. Then, a two- 
dimensional model of finite elements of plane deformation was developed to simulate the dynamic characteristics of 
cracking of the white layer. According to the hypothesis put forward, the contact area on the rail surface was plastically 
deformed under the interaction of the wheel and the rail; therefore, the energy-based J-integral criterion was introduced 
to evaluate the crack propagation behavior. 

Samples containing white layers on the surface of the railhead were cut from the rail, which used to be 
operated on the East Siberian Railway. The rail was made of pearlite steel in accordance with GOST51685-2013. This 
is grade 76KhfF steel: C 0.78%, Si 0.54% and Mn 0.8%, Cr 0.40%, V 0.035%. The mechanical properties of the rail 
were determined: 

- yield strength was 944 kN, 

- temporary resistance was 1287 kN, 

- elongation was 11.5%, 

- relative reduction was 31%, 

- impact strength was 18 J/cm’. 

The test temperature was +20°C. Brinell macrohardness on the rolling surface in places without white layers 
was 435, 485, 445, 465 and 494. This does not meet the requirements of GOST 51685-2013 (should be 352-405 
according to Brinell). Hardening of the railhead surface during operation is due to the combined action of a number of 
physical mechanisms [4—16]. For the calculation, the Young's modulus EF = 206 GPa was adopted, the shear modulus 
was 80 GPa, and the density was 7850 kg/m°* [3-10]. The standard metallographic procedure was performed after 
sectioning the samples in the direction parallel to the motion direction. After cutting, the sample was poured into resin, 
polished with silicon carbide sandpaper (grain size 1200—3000) and polished with 0.5 um microsilica. The cross section 
was observed under an optical microscope after etching with a 2% alcohol solution of nitric acid for 7 seconds. The 
mechanism of a white layer formation was investigated through the SYSWELD program, which uses an analytical model 
of the volumetric heat release of Goldak double ellipsoid. Three types of heat source are prescribed: 2D Gaussian, 
double ellipsoid and 3D conical Gaussian. After the geometric parameters of the heating model are determined and the 
maximum volumetric heat release value is set, an approximate calculation with constant thermophysical characteristics 
is performed. All material properties necessary for modeling are set in the form of piecewise linear functions. In the 
SYSWELD program, simultaneously with the heat problem, the metallurgical problem is solved. In the latter case, the 
calculations are based on the Leblond model, which describes the metallurgical transformation of one phase into 


another. The solution to the mechanics problems is reduced to the determination of thermal deformations. To do this, 
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you need to specify in the program the following: elasticity modulus, Poisson's ratio, thermal-expansion coefficient, as 
well as hardening curves for metallurgical phases. When cracking in a white layer under the moving contact conditions, 
observations were performed using ANSYS/LS-DYNA software, which is intended for finite element modeling. In 
studies on the dynamics of railway crews, the well-known algorithm of the named software FASTSIM with an elliptical 
contact area is widely used [25]. Compared with the geometric dimensions of the wheel and rail, the contact area is 
quite small (not more than 15 mm” [25]). Therefore, the model of wheel — rail contact can be simplified to the condition 
of deformation of the 2D plane [26, 27]. White layers in the form of an arc 2 mm long and 0.2 mm thick were 
constructed on the longitudinal section (the parameters were obtained as a result of the analysis of WEL studies on the 
rail surface [4—15]). In order to reduce computational costs and ensure accuracy in the model, the adaptive grid method 


was used. In the region of the white layer, very fine grids were used (Fig.1). 


FEM of white layer with cracks 
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Fig. 1. Finite element modeling (FEM) of white layer with three cracks (a) and transformation of element contour integral into 
surface contour integral: arbitrary contour around the crack tip (b); different contours of elements selected to represent independence 
from the calculation path of J-integral (c); basic diagram of 2D virtual crack closure method for calculating stress intensity factor (d) 


The established minimum element size was 0.02x0.02 mm’. The elements had a higher grid density with a set 
minimum element size of 0.010.01 mm’ at the tips and faces of cracks. This enabled to accurately capture the stress 
gradients. Outside the cracking zone, larger grids were used with a minimum element size of 0.1<0.1 mm’ in the wheel 
— rail contact surfaces. They gradually increased towards the far sections of the field. At first, the rail was combined into 
776,919 elements with 789,084 nodes, which changed during the simulation due to the use of the adaptive grid method. 
Such a model size is acceptable for obtaining accurate calculations of contact with rolling [19, 20, 28]. A constant 
vertical loading force of 13,000 N was applied to the wheel, which corresponded to the equivalent maximum Hertz 
pressure of ~ 1.2 GPa [25]. The wheel was set to a rotation speed of w 43.5 rad/s, which equals a travel speed of 72 
km/h with a friction coefficient of 0.3. It was studied how the properties of cracks in a white layer are affected by their 


length, changes in the angle of inclination, load pressure, and friction coefficient. 
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The crack growth that destroys the surface was analyzed under the assumption of linear elastic fracture 
mechanics [1—3, 17—19]. For a more accurate solution, we used the criterion of elastic-plastic fracture mechanics, the J- 
integral method, which is applicable for both linear and elastic plastic solutions [29]. It was introduced to study crack 
propagation behavior. For numerical calculations, the J-integral can be obtained from the solution in the far zone [29]. 
In the scope of work with finite elements, the element contour integral should be transformed into the surface contour 
integral, as shown in Fig. 15, where I is the contour of the curve that defines the boundary of the J-integral, and it is 
directed counterclockwise from the lower edge of the crack to the upper one. The J-integral can be additionally 


estimated as: 


du; 


J = f[(wdy -T, —ds), 
where w is strain energy density; 7; 1s the thrust vector; u; is the component of the displacement vector; ds is the length 
increment along the contour T° [29]. 

The J-integral evaluation is implemented in LS-PrePost as a post-processing tool. This application is also 
suitable for 2D modeling of plane deformation [30]. To calculate the J-integral around the crack tip, various contours of 
the elements were selected (see Fig. 1c). It is important to note that the unloading conditions for a path-independent J- 
integral are not proposed in accordance with the fracture mechanics [31]. The principle of expanding the J-integral 
methodology beyond the admissibility of linear elastic fracture mechanics was to idealize elastoplastic deformation as a 
nonlinear elastic deformation. The identity of the load characteristics — the deformation for elastoplastic and nonlinear 
elastic materials is known [30, 31]. The elastic-plastic material follows a linear unloading path with a slope equal to 
Young's modulus. Non-linear elastic material gets unloaded in the same way as when loading. Thus, the analysis 
suggests that the nonlinear-elastic behavior is valid for an elastoplastic material. Unloading occurs near the crack tip 
after the wheel has passed along the rail. The calculated values of the J-integral can be valid only until maximum is 
reached. Therefore, only the maximum values of the J-integral are discussed in the present paper. The stress intensity 
factors within the framework of the finite element scheme were calculated by the 2D method of virtual crack closure 
[29-32]. The forces y and x needed to connect the nodes c and d (see Fig. 1d) are denoted by F’, and T., respectively. 

It should be noted that the theory of linear elastic fracture mechanics is based on the assumption that there is no 
plastic deformation around the cracks. At the same time, according to [1-3], the use of the theory of linear elastic 
fracture mechanics is correct since the nonlinear deformation of the material is limited to a small area surrounding the 
crack tip. The reasoning that a plane problem can be considered in connection with a small contact area is incorrect. 
Under the plane deformation, this region is strip and infinite. And in this case, we are talking about an essentially three- 
dimensional problem. In [32], the complexity of the rail geometry and the boundedness of the region in which the main 
changes in the stress-strain state occur. The authors of this paper propose to consider a simplified axisymmetric problem 
for a multilayer coating in the form of a piecewise inhomogeneous layer. In [33], the stress-strain state of a multilayer 
coating was studied in the vicinity of the wheel — rail contact area. It was shown that the rigidity and thickness of the 
upper coating layer affect significantly the equivalent and contact stresses. At the fixed mechanical parameters of the 
coating, with an increase in the thickness of its upper layer, the values of maximum equivalent and contact stresses 
decrease. The authors of [34] applied thin multilayer coatings on the surface of a railway rail in the vicinity of the 
lateral wheel — rail contact. The stress-strain state of these coatings was investigated at various values of their geometric 
and mechanical parameters during rotary movements. 


Research Results. Fig. 2 shows optical images of the microstructure of the cross section of a white layer on 
the rail surface. 
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a) b) 
Fig. 2. Features of cracks in white layer on rail surface detected through optical microscopy: general view of white and deformed 
layer at 200-fold magnification (a); 1 - front edge crack, 2 - mid-position crack at 500x magnification (5) 


Fig. 2a shows a segment of a white layer on the surface. Its depth is 10 um. The microhardness value in the 
white layer is HV), = 989. This corresponds to the WEL features presented in [4-15]. The WEL thickness can reach 
200-300 um depending on the operating conditions [4—6]. The image shows a sharp boundary between the white layer 
and the deformed pearlite structure of the basis rail metal. The depth of the plastic deformation zone between WEL and 
the basic rail material is ~ 70 um. According to the estimates [5—10], on the pearlitic rails in operation, in the spot of the 
wheel — rail contact, plastic deformation has a high concentration creating a thin shear zone of 502,180 um. At the same 
time, hardness tests show that due to volume contact, the depth of plastic deformation can reach 1-10 mm [4—9, 27-30]. 
This causes an important evolution of the microstructure on a narrower scale, 1.e., to bend or rupture of cementite plates, 
to reduction of the inter-plate distance. The mechanical properties also change since the tensile strength increases and 
the fracture toughness decreases for cracks parallel to the aligned cementite plates. Two different cracks are located at 
the leading edge of the WEL (see Fig. 25). They can originate and propagate at the interface between the WEL and the 
deformed pearlite microstructure. A detailed study on other samples showed that an inclined crack appeared on the 
leading edge of the WEL following the flow direction of the rail material. Then the crack went down into the rail 
material along the layer boundary. Fig. 2b shows a crack propagating vertically in the mid-position within the WEL. In 
the samples studied, two different types of cracks are observed at the trailing edge of the WEL. They are obtained from 
different WEL sections, but are located on the same rail sample and are very close to each other. The first-type crack 
originates and propagates along the boundary between the white layer and the plastic deformation zone (see Fig. 25). 
The second-type crack crosses the interface with a deformed microstructure, but shows a tendency to curvature along 
the line in the direction of perlite shear strain. Such behavior of cracks is described in [4—7]. Indeed, the cracks studied 
by us are located in the same white layer; therefore, they were subjected to the same loads. However, their growth and 
development rates are different. This indicates a change in the state of stress in the leading, middle and rear position in 
the white layer. This means that the mechanism for the subsequent crack propagation (1.e., behind it) should be 
different. As shown in [35, 36], the plastically deformed structure of the pearlite region directly below the white layer 
can play a significant role in effecting the crack propagation. Cracks that are formed due to fatigue when in the rolling 
contact can be divided into two categories according to the place of their origination: beneath the surface and on the 
surface. 

Typically, cracks beneath the surface result from a strong vertical load in conjunction with material defects. 
Most cracks originate on the surface due to the wheel — rail interaction, as well as the transfer of a large load to a small 
contact area. The contact area is elliptical. It 1s relatively small, and at the same time, it supports the entire wheel load. 
Cracks, that arise as a result of fatigue when in the rolling contact and are the result of intense shear stress in the wheel 
— rail contact area, will grow when these stresses exceed the allowable tensile strength of the rail. It is also possible that 


cracks will advance to the top of the rails. According to [11, 12, 23], the heating rate of the surface layer of the rail, 
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which occurs during the passing of a train, can exceed 10°°C/s and reach the temperature at which austenite is formed. 
After passing the train under rapid cooling, austenite turns into martensite. The wheel — rail contact time is extremely 
short (milliseconds). In this case, the steel of the eutectoid composition is heated from room temperature to 727°C at a 
rate of 10°°C/s. 

Thus, it 1s practically impossible to measure actual temperature changes during the period under consideration. 
On the other hand, fast heating and hardening can be modeled and controlled in the SYSWELD program. For this, 
heating rates of 20—1000°C/s were used, which were obtained due to the parameter of the velocity of the heat source 
moving along the rail surface in the speed range 5—100 mm/s. 


Fig. 3 shows a fragment of modeling the mechanism of formation of a white layer in the SYSWELD program. 








a b 
Fig. 3. Simulation results of ae layer formation on the rail head surface in SYSWELD program: a fraction in the white 
layer, % (a); distribution of residual stresses in the white layer, MPa (b). 
According to Table 1, we can judge how the level of residual stresses depends on the parameters of the path 
and the hardening rate. 
Table 1 
Residual stresses value (compression and tension) along the path width 


and depth at a hardening speed of 15 mm 


Distance along the center of the path 0.65 | 0.77 | 086 122 | 154 | 25 3.0 
axis (depth), mm 





Centre Sent from path axis 0.14 | 165 | 1.92 | 219 | 259 | 3.16 | 3.91 5.2 
(width), mm 


1 


The white layer modeled in the program consists of martensite. The layer morphology is similar to that 
observed in the studied rail (see Fig. 2). The hardness of the simulated WEL is 670-810 AV. This is slightly lower than 
in rails (WEL ~ 725-1050 HV). However, current simulations represent only one phase transformation cycle. From one 
to five repeated heat treatment cycles were simulated on the surface of the rail head. It is found that hardness increases 
to 700-850 HV. Obviously, martensite with the same hardness as in the white layer of the rail can be obtained after 
several passes of the wheels. When modeling the formation of the white layer, it turned out that the temperature change 
during the wheel — rail contact and the distribution of residual stresses play an important role in this process. 


Fig. 4 shows J-integral maxima at the tips of the simulated crack. 
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Fig. 4. Maximum values of J-integral at crack tips in the leading, middle, and trailing edges of WEL 
(crack length for calculation is chosen as 0.01 mm; 0.03 mm; 0.05 mm; 0.1 mm and 0.2 mm, separately) 


It should be noted that the length of the crack is measured from the rail surface to the WEL depth. In the 
presented simulation version, cracks no. 1, 3 originate at an angle o and B (83.13°), the maximum load pressure is 1.2 
GPa. At this, all three cracks have a length of 10 um. At the top of crack no. 1, the value of the J-integral reaches a 
maximum of 0.321 N/mm. The maximum values of the J-integral for cracks nos. 2 and 3 are much smaller — 0.18 
N/mm and 0.088 N/mm, respectively. The analysis of Fig. 4 shows that crack no. 1 at the leading edge of WEL is likely 
to grow compared to cracks in the mid-position and at the trailing edge. This fact is in agreement with rolling fatigue 
experiments [10—18]. It is shown that the boundary (beginning) of the white layer on the rail surface and the basis metal 
has the lowest resistance to the formation of a contact-fatigue defect, while the fatigue life at the end of the WEL is 
approximately 3 times higher than in the center and at the starting point. Fig. 4 shows that the maxima of the J-integral 
around the crack tips show a tendency to decrease with increasing length or depth. This phenomenon is more significant 
for crack no. | than for nos. 2 and 3. As the crack length increases, the maximum value of the J-integral for all three 
cracks decreases, respectively, from 0.3214 N/mm to 0.222 N/mm; from 0.18 N/mm to 0.14 N/mm and from 0.088 
N/mm to 0.062 N/mm. This is consistent with the results of [19] since the deeper surfaces of the crack are less affected 
by normal and tangential contact loads. However, Fig. 4 shows that, for crack no. 2, when its length is 200 um and 
corresponds to the WEL thickness, the J-integral maximum tends to the opposite, which has never been reported. In this 
case, the J-integral maximum is the largest in comparison with other short cracks no. 2. It can be assumed that crack no. 
2 will rapidly propagate to the boundary between the WEL and the substrate. This is clearly seen in Fig. 25 (arrow 2). 
The significant difference is due to material irregularities between the perfectly elastic WEL and the elastic-plastic rail 
matrix. The calculation results show that for all three cracks the values of the J-integral decrease slightly when the crack 
length increases from 10 to 50 um. If the deepening of crack No. 2 is from 50 to 100 um, the change in the J-integral is 
relatively small. On the contrary, there is a significant increase in the maximum value of the J-integral if crack No. 2 
propagates through the entire layer to the boundary between WEL and the rail matrix. The calculation results show that 
for all three cracks, the values of the J-integral decrease slightly when the crack length increases from 10 to 50 um. If 
the deepening of crack no. 2 is from 50 to 100 um, the change in the J-integral 1s relatively small. On the contrary, there 
is a significant increase in the maximum value of the J-integral if crack no. 2 propagates through the entire layer to the 
boundary between the WEL and the rail matrix. 

The maximum value of the J-integral increases sharply with increasing load pressure from 0.8 to 1.5 GPa: 

- from 0.27 N/mm to 0.53 N/mm for crack no. 1; 

- from 0.17 N/mm to 0.47 N/mm for crack no. 2; 

- from 0.04 N/mm to 0.1 N/mm for crack no. 3. 
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This suggests that cracks are more likely to propagate at high loads. The maximum principal stress around the 
crack tip can be used as an indicator to assess the fracture properties and the crack path [19-26]. For example, if no. 1 
crack length is 10 um, the maximum principal stress at the crack tip is 1589 MPa. For cracks nos. 2 and 3 with the same 
length, different stress concentrations are obtained around the crack tips and the maximum principal stresses are about 
1506 MPa and 1261 MPa, respectively. An increase in length from 10 to 100 um causes a decrease in the maximum 
principal stress from 1506 MPa to 1206 MPa. When a crack propagates down to the boundary of the white layer, a 
sharp increase in stress concentration can be observed. This is clearly seen from modeling the distribution of residual 
stresses along the depth of the layer (see Fig. 3). So, according to J-integral calculations, under the same load 
conditions, cracks in the leading edge of the white layer propagate with greater probability than in the mid-position and 
at the trailing edge. This is consistent with crack growth. In this case, due to the inhomogeneity of the material, a 
complex stress field arises when the middle crack propagates at the interface between the white layer and the base 
metal. This causes a high concentration of stresses and thereby accelerates the crack propagation. It should be 
mentioned that this work did not consider the anisotropic characteristics of the material caused by severe plastic 
deformation in the transition zone. It was shown in [4—10] that this type of shear deformation causes elongation and 
alignment of cementite colonies of the surface layer of the rail head (see Fig. 2a). It 1s revealed that shear bands and 
regions with a high density of defects are found in the transition zone of the white layer, and some cementite plates in 
this region are destroyed [4-8]. Although linear elastic fracture mechanics does not reflect the detailed physics of the 
situation [1-3], it helps to understand the directions of the observed crack growth. To obtain more detailed information 
on the dynamics of crack formation in the white layer on the rail surface, stress intensity factors should be taken into 
account. At the same time, rolling contact cracks experience a mixed load corresponding to a sequence of tensile 
stresses followed by a shear cycle [4-20]. Therefore, stress intensity factors in the direction of maximum tangential 
stress can also be used to predict the direction of the crack growth. 

Discussion and Conclusions. It is established that cracks in the leading edge of the white layer (crack no. 1) 
propagate from the edge of the layer along the interface between the white layer and the basis metal of the rail. Cracks 
no. 2 (in the middle of the layer) grow vertically before reaching the plastically deformed matrix. Cracks no. 3 of the 
trailing edge show two developmental trends. In the first case, the crack extends along the interface; in the second case, 
it crosses the rail material and gradually aligns in the direction of deformation. 

Comparisons of the J-integral maxima have shown that under identical load conditions, crack no. 1 is most 
likely to have a high growth rate compared to cracks nos. 2 and 3. With an increase in the crack length, the maxima of 
the J-integral of all three cracks decreased. A great impact of the load pressure on the development of all the considered 


types of cracks is established. 
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Introduction. Loss of stability and buckling of a round plate may be observed if the plate is loaded on the lateral surface. 








Stability of a nonlinear elastic plate under lateral compression 
G. I. Volokitin 


Don State Technical University (Rostov-on-Don, Russian Federation) 


The solution to this problem is based on a bifurcation approach. In this case, a plate is considered as a nonlinear elastic 
body. It is important to choose the relation between stresses and deformations in sustainability problems of nonlinear 
elasticity. Simple laws of state (constitutive equations) were considered in early works devoted to this problem, for 
example, material of the “harmonic type” (Sensenig). 

Materials and Methods. Equations of neutral equilibrium for round plates made of Murnaghan and Blatz-Ko materials 
are obtained. Assuming a uniform initial deformation on the plate, the stability problem is considered. Strict three- 
dimensional neutral equilibrium equations provide exploring related forms of equilibrium taking into account physical 
and geometric nonlinearity. Derivation of these equations is based on the application of the theory of superposition of 
small deformation on the final one. 

Results. Progress in solution to the corresponding secular equation (with non-linear parameter entry) for practically 
important laws of elasticity of Murnaghan and Blatz-Ko is possible using the numerical methods only. The developed 
method for calculating bifurcation values of loading parameters makes it possible to analyze the effect of nonlinearity. 
Discussion and Conclusions. The influence of physical and geometric nonlinearity on the upper critical value of the 
initial deformation parameter is explored. The results obtained can be used under the assessment of reliability of elastic 
third-order moduli for various physical materials. Data on these moduli is still scarce. The numerical research has 
shown that the constants given in some sources should be treated with caution. The use of elasticity moduli in the law of 
state of Blatz-Ko is also discussed. 


Keywords: finite deformation, stresses, buckling, plate, stability. 
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Introduction. Currently, the study of new relatively simple and adequate laws of state for various materials 
that require considering nonlinear deformations is an actively developing area of continuum mechanics [1—10]. In the 
framework of the theory of superposition of small deformation on the finite one, three-dimensional equations of neutral 
equilibrium are derived for the materials of Murnaghan and Blatz-Ko. Based on these equations, an example of end 
buckling of a uniformly compressed round plate is considered [11, 12]. 

Materials and Methods 

Equations of neutral equilibrium. Let 7, @, z be the cylindrical coordinates of a point in an undeformed 
state. We assume that the initial deformation of the body is determined by the radius vector R: 

R= R(r)e, + dzi, , (1) 


where R(r)=ar, parameters a, d are constants, e , e, 1, are basis vectors. 


For the coordinates of this point in the initial deformed state, we have: 
R=ar, D=o0, Z=dz. 
Therefore, the strain gradient, the Finger strain measure, and its principal invariants are determined through the 
relations: 
VR=VR’ =a(e,e, + Cpe) + dizi,, F= VR’, 
I,=2a°+d*, I,=a'+2a’d’, I,=a‘d’. (2) 
From now on, V, V_ is the nabla operator in the metric of the undeformed and initial-deformed state: 


V=VR7-V. 
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We will use the neutral equilibrium equations proposed by A. I. Lurie [13]: V@=0, where the tensor © is a 
linear differential operator over the vector of additional displacements W. The expressions of the components of this 
tensor concretized with account of the laws of state of Murnaghan and Blatz-Ko were obtained in [14, 15]. In the 
representations of the tensor ©, the components are some functions defined as a result of solving the boundary value 
problem of the initial deformation. 

Research Results. To consider the bending form of the plate equilibrium bifurcation, similarly to [11, 12], we 
accept the following additional displacement vector 

W =Uu(r,z)e. + W(T, Z)i; . (3) 

Considering (1), (2) and (3), Lurie tensor has the form: 

© = Ae_e, + Be,e, + Ci,i, + Gi,e, + He, i, . 


Here, 
A=A, HA p= B 0 eae as ,~C=G “40,4, 
Or Oz * Or Oz Or OZ 
G= 3 7. He ee ee 
OZ Or OZ Or 


We note that the structure of these operators is typical in the stability problems of cylindrical nonlinear elastic 


bodies. Omitting the expressions of the remaining coefficients, we give, for example, formulas for A,. In the case of 


Murnaghan law of the state [13, 14] A, is expressed through the relation 


2 
Henceforward, 4 and pt are Lame elasticity moduli, v,, v,, v, are the third-order elasticity constants. For 


Blatz-Ko material [15]: 


A 
=9{1-A y that a(y, 420%". 
d 2 


A= 2u(1—B) 7 Pl gor aE B Fo) 


1 Zi Ji, 
where f is the refining elasticity modulus. 
The equation of neutral equilibrium is equivalent to the system of differential equations with respect to the 


2 2 2 
A, aa [4+4 +—G ; a Foamy 
Or’ ror r d oroz d 


2 Z 
(A, + H, )— g (4 4),[ 2 PAB), ec on 2G 


components of the vector W : 





(4) 





dz r or ror) d > @& 


Assuming that the plate is loaded with uniform pressure along the lateral surface, we supplement the 
differential equations (4) with the equilibrium conditions at the boundary: 


u(r,z)|,_, ie YZ |r=r, = 0. 





At the ends, 1.e., at 
z= 45 © =0, 0, =0. (5) 


We apply the substitution 


‘ =X, (Zz) (K,7); (6) 


w=Z_(z)J,(K,r). 
Here, n=1,2,..., J,, J, are Bessel functions of zero order and first kind, and numbers k,r, are null functions J,(x). 


Let us assume x =k,r and use the equations 


2 
OE 2 gS), AED 2 prs, Po 2420) = 4,00, 
dr dr dr 
and well-known identities for Bessel functions 
1 | | 


ee, 8) =, 
X X X 


Having completed the above steps, we obtain that the variables in the differential equations (4) and under the 
boundary conditions (5) are separated. We arrive at the following boundary value problem for ordinary differential 
equations: 
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A G.-C, 


X"-k° > X -k, -~—Z' =0, 
2 G, 
G G,+C “) 
Z" -~k° “7 +k, -~— X' =0. 
C; C; 


k,C,X,, (£h/2)+C,Z} (£h/2) = 0, 

k,G,Z, (£h/2)-G,X) (£h/2) =0. 

Two homogeneous linear ordinary differential equations (7) and four boundary conditions (8) result in the 
eigenvalue problem with a nonlinear occurrence of the parameter. In this problem, such a parameter is the 


(8) 


quantity (1—a) . The system (7) is ended in the standard form: 





Vi =)>, 
A G,+ 
’ 2 2 2 2 
V> =k; —— ark 
Al, G, 
V3 =Va> 
G G,+C 
’ 2 1 2 2 
¥, =k, 3 -K, J: 
C; C; 


The following notations are accepted here: argumentt=z, (¥333)33),) =(X,;X1;Z,;Z/) . Let the 

fundamental system of solutions be the following four vectors y, =(¥1,3 ¥2;3 ¥33 Va i . Then, the general solution to the 
4 

system is given by (Xx LA, y = »<¢,Y,; » where ¢, are arbitrary constants. For example, we set the initial data 
‘ js] 


for z=—h/2 by the columns of the fourth-order unity matrix. Then we solve numerically the Cauchy problem with 
these initial conditions. As a result, on the right-hand side, we obtain the values of the basis functions, 1.e., 
vectors y,(+//2) . Using the boundary conditions (8), we arrive at a homogeneous system of linear algebraic equations 
with respect to ¢,, ¢,, ¢3,5 ¢,: 

4 

4,6, =0, where i= 1, 2, 3, 4. 

j=l 


Here, the coefficients are the elements of matrix A: 
a; =k, GY; (—h/2)+ Cra; (—h/2) » A, =K, CN; (1/2)+ Cra; (h/2) , 
a,, =k,G,y;,(-h/2)-G,y,,(-A/2), a4, =k,G,y;, (4/2) - Gy, (4/2). 
The homogeneous system has a nontrivial solution if the condition is met 
detA =0. (9) 
The determinant expression includes the load parameters a, d , Bessel null-functions, as well as the elasticity 


moduli 4, w, v,, v,, Vv; (for the Murnaghan material) or 1, “, 6 (for the Blatz-Ko material). 


The parameters a, d , which set the initial deformation, are interconnected. The axial force acting on the cross- 
Me 

sectional area is determined by the relation [13]: O =22[|0,RdR, where oa, is the physical component of the stress 
0 


tensor. Since the initial deformation is assumed in the form (1), the Cauchy stress tensor T and the Finger strain 
measure F are coaxial. 
Moreover, the stress tensor 1s constant: 


2 - 7 ns 
T = ae (a°(e,e, +e,e,)+d’i,i,)—c" (a*(e,e, +e,e,) + di,i,) +c (ee, +e,€, + i,i;)) 


The ends of the plate are free of load, therefore o, = 0. This implies the condition relating the coefficients of 


the Finger law 
oD 
c —cd’ + =0. (10) 
d 
Considering the law of state, we can write a specific expression of the condition (10), which establishes a 
connection between a and d. So, for the Murnaghan material, we get the condition: 


3 3 
(4 ay, Ja {2424 Ah-sy, —4yv, +(y, 1 2y)a" Ja — 
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ov, 9 3 15 
“3h— 2+ +2, + (2h 3y —4v, — a? Hy, $Y a’ =0. 
In particular, if v, =v, =v, =0, then we get 
fon (1-a’). 
A+2u 


When considering the Blatz-Ko law of state, d is expressed by the formula 


d=a, 
Thus, through setting the elasticity moduli and Bessel null-functions, we find the bifurcation values of the 





initial strain parameter a, from (9). 
Note that the value A =1-—a, is a small parameter for relatively thin disks. So, in the classical theory of plate 
buckling, the critical value of a, is determined by the formula 
> 2 
_|_ (3,8317)° | A 
120 +v) lr, J 


where v is the Poisson's ratio, and 3.8317 is the first root of the Bessel function J,(x) . 


* 


When solving the initial boundary value problem, it is assumed that the Signorini’s perturbation method can be 
applied. At this, we assume that the coefficients of the operator © depend on a small parameter A in a power-law 
manner. This means that the boundary conditions for the incremental displacement components A are specified by the 
matrix A = A(A). In this problem, a partial eigenvalue problem — the determination of the lowest eigenvalues — has 
physical meaning. The higher degrees A , which contain the matrix elements, slightly affect the values of the smallest 
roots of the secular equation (9). If we restrict ourselves to the linear theory under solving the initial problem, then 
A(A) is a regular binomial [16]. 

Numerical experiments under studying the stability of nonlinear elastic bodies of not too large relative 
thickness confirm this conclusion. Therefore, a characteristic equation with non-linear occurrence of parameter (9) can 
be replaced by the characteristic equation of the linear operator. Iterative processes that converge to one eigenvalue, 
where a number close to the value in the theory of plates is chosen as the null approximation, can be applied. 

Discussion and Conclusions. Using the equations obtained above, a numerical analysis of the influence of 
physical and geometric nonlinearity on the value of the upper critical parameter a, 1s performed. The calculations are 
implemented in the Matlab environment for various options of specifying elasticity moduli, relative plate thickness, and 
waveformation number 7 . It is found that in all cases of loss of disk stability, the first-mode buckling that corresponds 
to the minimum critical value of the parameter, which corresponds to the root of the Bessel function 3.8317 [17], 


answers the smallest critical parameter value a, . 
Table 1 shows the critical parameter values a, for plates with a relative thickness of 0.05; 0.1; 0.15; 0.2; 0.25; 


and 0.3. 
Table 1 


Critical parameter values a, for plates of various relative thicknesses 


h/r Relative plate thickness 


| 
0.9972 0.9867 0.9735 0.9549 0.9300 0.8914 
0.9967 0.9869 0.9690 0.9410 — — 
0.9965 0.9856 0.9665 0.9372 — — 


0.9881 — — 0.9965 0.9613 0.9861 
0.9929 0.9729 0.9341 0.9114 — — 

| | toms | os | ton | osms | toes | oss 
0.9985 0.994] 0.9871 0.9775 0.9651 0.9539 
0.9985 0.994] 0.9863 0.9749 0.9594 0.9396 
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The following numbers indicate the materials listed below. 

1. Steel Rex 535 (4 =1.09, w=0.818, v, =-1.75, v, =-2.40, v, =—-1.69). 

2. Steel SOHGSM2F (1 =1.129, w=0.803, v, =—-2.8, v, =-2.7, v, =—1.87). 

3. Steel Hecla37(A=1.11, w=0.821L v, =—-3.58, v, =-2.82, v,=-1.77). 

4. Steel Hecla ATV (A=0.87, w=0.716, v, =0.34, v, =—-5.52, v; =—-1.0). 

5. Beryllium bronze (A=/.042, u=0.49, v,=—4.0, v,=—-1.7, v, =—-0.6). 

6. Organic glass (A=0.39, u=0.186, v, =—0.078, v, =—0.07, v; =0.047 ). 

In the first six cases, the Murnaghan law of state is considered. The last option presents the results for the 


material of Blatz-Ko. The top number in the table cell refers to the case in which physical non-linearity is not taken into 
account, Le., v, =v, =v; =0 in the Murnaghan law; £6 =0 in the Blatz-Ko law. The second (lower) number takes into 


account physical nonlinearity. In the Murnaghan law, data were used for the third-order elasticity moduli from [14] in 


units 10” Oe 


cm 
In the last row of Table 1, there are the results for the Blatz-Ko law, which was chosen in a simplified version 
(the Knowles-Sternberg equation): the Poisson's ratio was taken equal to 0.25, and the refining module 6 =0.45. A 


dash means no critical value a, was found. 


The analysis of the results provides drawing some conclusions. At small relative thicknesses of the disk, the 
exact theory and the linear theory of plates give the same critical parameter values a,. Geometric nonlinearity has 


marked impact at relative thicknesses greater than 0.1. Physical nonlinearity is even more pronounced. However, care 
should be exercised in choosing third-order elasticity moduli. For example, in the fourth version (Hecla ATV steel) and 
for thin plates, no reliable critical values of the parameter a, were found; although acceptable values indicating a loss of 


stability are observed for the same Lame elasticity moduli. As in the stability problem for a nonlinear elastic sphere 
made of Blatz-Ko material [15], the constant 8 weakly affects the critical value a, . 
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Development of design methodology of technological process of ball-rod 
hardening with account for formation of compressive residual stresses 





M. A. Tamarkin, E. E. Tishchenko, S. A. Novokreshchenovy, S. A. Morozov 


Don State Technical University (Rostov-on-Don, Russian Federation) 


Introduction. The study results on the multicontact shock-vibrating machining using a ball-rod hardener are presented. 
Methods of surface plastic deformation processing, their advantages are described. The tool circuit diagram is given. 
Features, technological advantages, and the application area of the ball-rod hardening are specified. 

Materials and Methods. When conducting theoretical studies on the processing, factors, which affect the quality of the 
surface layer of the machined parts, were established. Dependences are given for calculating the surface roughness, the 
hardened layer depth, and the deformation ratio under ball-rod hardening. While studying the generation of residual 
stresses, the dependence for calculating the residual stresses generated in the surface layer of the machined part was 
specified. 

Results. The experimental findings of the processing requisite for verification of the adequacy of the theoretical models, 
as well as the routine of the experiments, are presented. A table and graphs clearly confirming good convergence of the 
theoretical and experimental data are given (the difference does not exceed 20%). Residual stresses in the surface layer 
are compressive which enables to predict high performance properties of the machined parts. The value of residual 
stresses on the workpiece surface is in the range of 130 + 200 MPa. The depth of compressive residual stresses is in the 
range of 0.9—1 mm. The fatigue characteristic variation, the ultimate stresses of the cycle in depth, which affects the 
endurance limit, is calculated. It has been established that the processing of workpieces by a ball-rod hardener provides 
increasing the ultimate cycle stress under repeated loading by 27—35%. 

Discussions and Conclusions. The design methodology of technological process of ball-rod hardening can be used 
under the development of production at the machine-building enterprises. In accordance with the recommendations, the 
limits of the required quality parameters of the workpiece surface layer are set; the parameters of the ball-rod hardener, 
the interference fit and the radius of rod sharpening are selected. Quality parameters of the surface layer are calculated. 
Correction of the selected modes and re-calculation of the parameters of the machined surface are carried out until all 
the specified characteristics are located within the required limits. 


Keywords: ball-rod hardening, surface roughness, hardened layer depth, deformation ratio, residual stresses. 
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Introduction. The quality of parts depends on many different parameters, one way or another manifested in 
the process of their manufacture and affecting their life cycle. A special role in the formation of surface quality is 
played by finishing processing methods, the use of which provide improving the basic operational properties of machine 
parts, such as contact stiffness, wear resistance, fatigue strength and durability, etc. 

Surface plastic deformation (SPD) processing methods allow solving the problems of increasing the product 
life cycle for parts from various materials of any configuration. In addition, such methods provide localized processing, 
which can significantly reduce production costs. The characteristic surfaces of the parts processed through the local 
SPD hardening are stress concentration zones (holes, slots, bevels, sampling, screwthreads, fillets, grooves, etc.); 
unreinforced areas of the surface of parts that have undergone general hardening in vibration, shotblasting and others 
installations (under clamps, in pockets, holes and other difficult-to-access areas for the processing environment); places 
for mechanical refinement of parts, etc. 
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Methods of processing SPD, including local ones, can significantly reduce the surface roughness of the 
workpiece and increase its physical and mechanical properties. In addition, the SPD methods make it possible to 
remove residual stresses generated by various defects in a part, which have a significant effect on the 
physicomechanical properties of the material from which the part is made. It is known that tensile residual stresses can 
significantly reduce the strength characteristics of a part, and vice versa, compressive ones can improve them. 

To carry out the process of local processing of SPD parts of simple and complex configurations, including 
those having a small vertical difference in height, a special device, a ball-rod hardener (BRH), was invented at the 
Engineering Technology Department, Don State Technical University, under the guidance of Prof. Babichev A.P. BRH 
is a multi-contact vibro-impact tool, the processing of which is based on the surface-plastic deformation. This 
processing method has several advantages, such as the ability to handle stress concentrators of large-sized products, 
products of irregular shape, low-rigid products. The method combines the technological capabilities of vibratory 
finishing and hardening processing and embossing. 

A pneumatic hammer, on which the hardener body is fixed, is used as a power drive. The striker of the power 
drive 1 hammers with a frequency of about 42 Hz on several layers of steel balls 4, which are transmitted to the set of 
round rods 2 installed in the collet clamp 6 (Fig. 1). Several layers of balls allow rods with a spherical sharpening of the 
working surface to copy the shaped profile of the workpiece 5. 
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Compressed 
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Fig. 1. Diagram of multi-contact vibro-impact tool (BRH): 
1 - power drive, 2 - set of round rods, 3 - hardener body, 
4 - steel balls, 5 - workpiece, 6 - collet clamp, 7 - elastic element 


Main Part 
The processing method under consideration is studied by a fairly large number of authors [1—10]. When 


studying the technological capabilities of the BRH processing, it was found that the impact energy of the drive, the 
diameter of the sharpening of the rods, the number of rods in the set, and the advance of the device along the working 
surface, have the greatest impact on the quality of the workpiece surface layer [1, 2]. In the work of Shcherba L.M. [1], 
a methodology for designing the process of treatment by a ball-rod hardener is developed. 

Dependences are presented for determining the quality parameters of the surface layer: roughness of the treated 
surface, the degree of deformation, and the depth of the hardened layer. In the work of Isaev A.G. [2, 5], these 
dependences are refined and described below. 


The surface roughness is determined by the formula 


Re = 0.03, {2 1 
TV DN. HBP? ’ ) 


where £', is indenter impact energy; N is the number of rods in the set, HB is hardness of the work material 


according to Brinell; n is the device efficiency depending on a number of factors (interference during processing, the 
number of layers of balls, etc.), D is the diameter of the spherical sharpening of the hardener rod (indenter). 


The deformation ratio is presented as: 


e=1137 sl (2) 
D>. N- HB” 
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The depth of the hardened layer is determined from the following relationship: 








(3) 





In the study of Yu.R. Kopylov [3], the dependence of the determination of residual stresses for vibration 
finishing and hardening treatment is given. It seems possible to apply a similar technique to determine the residual 


stresses under the BRH processing through introducing the coefficient £, that considers the features of the generation of 
internal stresses under the BRH processing. The coefficient k, is determined during experimental studies. After 


transformations, the dependence will take the following form: 


5 


E l-v; l-v: (5) 
ana | ae 
7) ul 


2(D+d) 








: (4) 





where k, is the coefficient that considers features of the generation of internal stresses during the BRH processing 
(determined experimentally); k, is the coefficient that takes into account the ability of the material to hardening; £,, is 
indenter impact energy; D and d are the diameters of the rod sharpening and the indent on the workpiece, 
respectively; vy, and v,, are Poisson's ratios for the workpiece and the rod; Ey and E,, are elasticity moduli for the 
workpiece and the rod, respectively. 

To verify the adequacy of dependence (4), experimental studies on the generation of residual stresses under the 
BRH were carried out. Residual stresses were measured according to the standard method used by the factory 
laboratory, Rostvertol, Russian Helicopters JSC, using an automated test bench for residual stresses ASCON-3-KI on 
samples of B95 aluminum alloy. The rectangular samples 20010020 mm were processed under various modes with a 
tightness of 1.5 mm and 4.5 mm. Sets of rods with a radius of spherical sharpening of 4 mm and 8 mm were used. The 
residual stresses were determined according to the Davidenkov method through etching surface layers from plate 
samples cut from rectangular work samples. The measured values were compared to the theoretical values calculated 
from the formula (4). The calculation results are presented in Table 1. 

The coefficient value k, is determined from the results of preliminary experimental studies through comparing 
the magnitude of the residual stresses calculated by the dependence (4) and obtained experimentally. It varies in the 


range of 1.5—1.7. For calculations, the value of 1.6 is taken 
Table 1 


Comparison of the results of theoretical and experimental studies 


Initial data Difference of theoretical 
Theoretical values, , 
Rod sharpening Experimental data, MPa data from experimental data, 
Interference 
radius 
4 





Fig. 2 and 3 show graphs of the dependences of the residual stresses on the interference under processing. The 
solid line shows the theoretical curve. The squares indicate the results of experimental studies. 
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Fig. 2. Dependence of residual stresses on interference under processing. 
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Fig. 3. Dependence of residual stresses on interference under processing. 
Spherical sharpening radius of the rod R = 8 mm 


For the values of residual stresses obtained as a result of studies, the changes in the fatigue characteristic, the 
ultimate stresses of the cycle in depth, were calculated using the formula given in [1]. It is known that the endurance 
limit value is affected by the value of the average stress of the cycle, which largely depends on the level of residual 


stresses of the surface layer. 


The graphs are presented in Fig. 4—7. The solid line shows the ultimate stresses of the cycle, the dotted line 


shows the distribution of residual stresses. 
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Fig. 6. Change in limiting stresses value of the cycle along the Fig. 7. Change in limiting stresses value of the cycle along the 
workpiece depth after BRH processing: R = 8 mm, interference workpiece depth after BRH processing: R = 8 mm, interference 

is 1.5mm is 4.5mm 
Conclusion 


Based on the results of the research, the following conclusions can be drawn. 

1. The residual stresses in the surface layer are compressive, which enables to predict high performance 
properties of the finished workpieces. 

2. The magnitude of the residual stresses on the workpiece surface is in the range of 130—200 MPa, which, 
according to the data of Rostvertol PJSC, approximately corresponds to the magnitude of the residual stresses after the 
vibration hardening treatment. 

3. The depth of the compressive residual stresses is in the range of 0.9—1 mm, which, according to the data of 
Rostvertol PJSC, is approximately two times greater than during the vibration hardening treatment. 

4. BRH processing of workpieces provides increasing the ultimate stress of the cycle during cyclic loading of 
the workpiece by 27-35%. 

5. Processing modes have a slight impact on the magnitude and depth of the residual stresses. 

6. It is established that the difference between the results of theoretical and experimental studies does not 
exceed 20%, which indicates the adequacy of the obtained theoretical dependence for calculating the residual stresses 
under processing by a ball-rod hardener. 

Studies on residual stresses with BRH are presented in the works of A. Shvedova [4, 6—9]. In [4], computer 
simulation of the process was performed using the Ansys software package. The obtained values of residual stresses 
correlate with the results of the above experiments. 

Based on the research results, the following methodology for designing a ball-rod hardener treatment is 
proposed. 

1. The limits of the required quality parameters of the workpiece surface layer are set. 

2. A pneumatic hammer is selected. Then, a nozzle is selected considering those at the tooling engineer's 
disposal. For processing small areas, nozzles with a small number of rods are recommended; for processing areas of a 
large area — with a large number of rods. The number of layers of balls is selected depending on the height of the drops 
or the radius of curvature of the work surface. 

3. Interference during processing and the radius of sharpening of the rod are selected. For harder materials, 
higher machining interferences and smaller sharpening radii are selected. 

4. We assign a processing time of 10-15 seconds for the area of the bundle of rods. Considering the feed value 
recommended above, the number of tool passes along the workpiece surface is selected, and in most cases it is desirable 
to use single-pass machining. If there are several options for combinations of processing modes that enable to obtain the 
specified characteristics of the surface layer, the one, whose total processing time for a particular workpiece is shorter, 
is selected. 

5. Then, the arithmetic average deviation of the surface roughness, the depth of the hardened layer, the 


deformation ratio, and the residual stresses are calculated according to formulas 1-4. 
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6. Based on the calculation results, the selected processing modes are adjusted. Then, the parameters of the 
machined surface are calculated again, and so on, until all the specified characteristics are located within the necessary 
limits. 

The above technique can be used under the production design of the BRH treatment. 
The above engineering recommendations were used in the implementation of the BRH treatment at 
Rostvertol PJSC. 
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Introduction.The paper is devoted to the comparison and evaluation of the quality of a car using the method of mathe- 
matical simulation of the integral performance index. The work objectives were to develop a mathematical modeling 
technique for the integral performance index, to draw up a step sequence of quality assessment, to analyze most com- 
mon brands of passenger cars on the domestic market, to sum up and validate the result obtained. 

Materials and Methods. The sequence of simulating the integral performance index is proposed. Numerical modeling of 
the integral performance index is carried out on the example of six most common brands of passenger cars on the do- 
mestic market. 

Results. A method of modeling the integral performance index is developed. A step sequence of the analysis is de- 
scribed. An additive form of determining the integral performance index which provides combining single, unrelated 
quality indices into a complex indicator in the process of comparison 1s proposed. 

Discussion and Conclusions. The proposals presented can be used in the process of assessing and diagnosing the com- 


petitiveness of not only cars, but also of other products with a large set of independent quality indicators. 


Keywords: modeling, integral index, quality requirements, motor car, performance specifications, additive method, 


weight coefficient. 
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Introduction. To obtain, systematize, describe and apply accumulated knowledge and information in any 
field of activity, a universal simulation technique is most widely used. Mathematical modeling is systematically sup- 
ported and updated in many fields of science. Mathematical modeling provides combining accumulated knowledge 
and real processes including mental ones [1], combining the ability to reproduce the properties of a real or created 
object, process or phenomenon using another object, process or phenomenon. 

Quality is a combination of consumer properties unrelated to each other through technical characteristics that 
must satisfy current and future needs [2]. At the same time, a close relationship between quality and requirements can 
be traced. Requirements are specific features and conditions that correspond to consumer preferences throughout the 
life of the product. Product requirements are laid at the design stage [3, 4]. 

However, quality and requirements have some inconsistency [5, 6]. This is most pronounced in the discrep- 
ancy between the declared technical qualities of the goods and the requirements of consumers. This is due to the con- 
stant change in consumer requirements depending on technology development, as well as the financial and cultural 
condition of the population. 

Depending on the number of characterized properties and on their influence on the quality of the product, qual- 
ity indices are divided into simple and complex. A simple index enables to characterize one property or one dependence 


of quality on the technical and operational characteristics, while a complex index combines several characteristics. 
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There is a concept of direct and inverse indices. An increase in direct indices causes an increase in quality; an increase 
in inverse indices causes a decrease [7]. 

After analyzing the classification of product quality indices (Fig. 1), we can conclude that the integrated in- 
dex is a special case of a complex index. An integrated quality index characterizes the ratio of the overall useful effect 


during the operation of a given product to the cumulative costs of its purchase, operation or consumption [8]. 


Product quality indices 


[Simple J} 1 {copes 


Fig. 1. Classification of product quality indices 









Relative Integrated 


Simulation technique for the integrated index of product quality. 

The definition of a quality index by an example of a car is carried out in the sequence: 

¢ goal setting; 

¢ selection of prototypes for comparison and choice of the optimum alternative; 

¢ selection of parameters characterizing the product from the standpoint of the consumer; 

¢ determination of averages of quality indices; 

¢ combining averages into one index for each brand; 

¢ analysis of the integrated indices obtained for all prototypes; 

¢ decision-making on product quality control. 

The greatest difficulties are the selection of characteristics and obtaining the averages of these data [9] since 
there are no common techniques for forming a list of quality indices and transforming each index into a numerical 
form. To determine the integrated quality index, its additive form is widely used, 1.e., the weighted average summa- 
tion, which allows combining simple, unrelated quality indices into a complex index. It is used in cases when deci- 
sions are made on the economically best option if several prototypes are available on the market. That is, this tech- 
nique is most suitable for modeling an integrated index of car quality. 

The formula for calculating the complex integrated index [10]: 

K; = ie WG * Xi 
where a, is weight coefficient of the i-th parameter; X; is the averaged quality index for the i-th parameter; n is the 
number of parameters to compare. 


To determine the weight coefficient, an expert method is used according to the expression: 
a, =1.0. 
i=l 


The algorithm of the calculations used can be traced on the example of choosing the best option for buying a 


car (Fig. 2). 
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Fig. 2. Popular sedans at the domestic market 


The analysis was carried out on the bestselling cars at the domestic market. The hierarchy is implemented by 
price characteristics. 

Study Results. For the purity of the experiment, six sedan cars were compared in the same price category. 
There are several tens of quality indices; in the calculations, an integrated index was calculated for 6 groups of basic 
(from the consumer's point of view) quality indices. The initial data for the calculation are given in Table 1. In prac- 
tice, the use of such quality assessments can take into account a different number of characteristics, even considering 
the fact that many of them have a rather contradictory effect on quality and, in their physical nature, are very difficult 


to quantify. 
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Table | 
Source Data 
Model 
Car parameter Lada Hyundai Skoda Kia Lada Vesta Volkswagen 
Granta Solaris Rapid Rio Polo 
Price, thousand Rub 580 750 680 640 600 690 
Mileage, thousand km 50 q 5/ 26 a3 50 
Fuel rate, | Tez 8.3 8.5 6.4 8.7 8.0 
Power, kWt 72.07 90.50 80.91 73.50 78.00 80.91 
Clearance, mm 170 160 136 160 171 163 
Acceleration time, s 11.5 11.0 10.0 13.1 12.8 11.7 








An intermediate stage of the calculation is to obtain the averages of the 6 considered quality indices for each 


model. The weight coefficient is determined by the expert evaluation method, where the human factor (personal prefer- 
ences) also plays an important role in the assessment and significantly affects the final value of the integrated quality 


index of the sedans compared. The calculation results are presented in Table 2. 












































Table 2 
Calculation Results 
Model 
Quality index ; 
iad Geant Hyundai Skoda Kia Lada Volkswagen 
Solaris Rapid Rio Vesta Polo 
Rank by sales l 4 6 3 pi 5 
Weight coefficient 0.11 0.15 0.25 0.15 0.14 0.20 
Int ted lit 
ee 0.32 0.32 0.71 0.48 0.33 0.49 
index 
Rank by quality index 5 5 l 3 4 2 
Marketability 55-10 42:10° 10° 75:10° 55-10 71:10° 
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Thus, the maximum quality integrated index belongs to Skoda Rapid and it is 0.71. At the same time, the val- 
ues of the quality factor of the other compared cars that are in the same price range differ slightly. The values of the 
competitiveness indicator of the compared sedans are obtained by the ratio of the quality index and the cost of the car. 
This is also maximum estimator for Skoda Rapid, which is explained by the fact that the buyer of this car for one mone- 
tary unit gains more quality than buyers of the prototypes compared. 

Discussion and Conclusions. Recording information about an object of research or design requires using 
mathematical simulation to store and transmit it in space or in time. The simulation is aimed at building, improving, 
studying and applying models of actually existing or designed objects for the subsequent assessment of their competi- 
tiveness. The integrated quality index is one of the unified indicators that can implement a comparison of indices even 


with a conflicting effect on quality. 
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Mobile machine design through dynamic load simulation on their drive units 
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Introduction. The mobile machine design is impossible without considering the vibration parameters of their units. This 
requires the development of specialized dynamic models that take into account the probabilistic nature of these 
parameters. The root cause for the occurrence of vibration effects is the profile irregularity of the mobile machine path, 
and the variability of physicomechanical characteristics of the soil. Problems that consider these features are solved 
linearly with sufficient accuracy; but in multidimensional dynamical systems, such an approach is unacceptable due to 
the presence of a large number of interrelationships. The work objective is to conduct a comparative analysis of the 
efficiency of existing calculation methods of statistical characteristics of uncorrelated external actions as applied to a 
mobile machine presented as a multidimensional dynamic system with actions having different correlations. 

Materials and Methods. External actions in the multidimensional dynamical systems are considered in a matrix form. 
When calculating statistical characteristics, intercouplings in the spectral density matrices are taken into account. The 
elements of the main and secondary diagonals are determined; the correlations between the effects are taken into 
account. These features significantly complicate the calculations. So, to get matrices of uncorrelated actions, the matrix 
of external actions is reduced to a diagonal form. 

Results. A numerical comparison of spectral densities and intensity of the mobile machine oscillations under variation 
of speeds and nature of the soil fertility microprofile was carried out using various methods of calculation. Certain 
characteristics of spectral densities and oscillations of mobile machines of agroindustrial complex enabled to develop 
recommendations on the practical application of the presented dependences for designing this machinery. 

Discussion and Conclusions. The results of solving the matrix of spectral densities of external actions by various 
methods are presented on the diagram of spectral oscillation velocities. The analysis of characteristic curves has shown 
that the identical results, regardless of the calculation method, are obtained only for machines with weak functional 
relations under the uncorrelated external action of the soil fertility. For some cases, the resonance machine speeds are 
set. The effect of irregularities of the soil fertility on the oscillation intensity of the machine units and the dispersion of 


loads on the units is shown in the graphical representation. 


Keywords: multidimensional dynamic system, spectral density matrix, uncorrelated external influences, inverse matrix, 


attached matrix. 
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Introduction. When designing reliable mobile equipment, both static dynamic loads should be considered. 
Vibrations of the machine and its components affect significantly the characteristics of dynamic loads. Accounting and 


determination of vibration parameters require the development of specialized dynamic models that consider these 
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parameters as regards to probabilistic analysis. The root reason for the vibration effects is due to the uneven profile of 
the mobile machine path and the inconsistency of the physical and mechanical characteristics of the agricultural 
background. Problems that consider these features are solved linearly with sufficient accuracy, which is unacceptable 
when solving multidimensional dynamical systems with a large number of interconnections. 

Materials and Methods. Quality of the technological process performed by a mobile agricultural machine 1s 
affected by both the characteristics of its working bodies [1—5] and the characteristics of the undercarriage system [6, 
7], aS well as the microrelief of the agricultural background [8, 9]. In general, external actions on a machine are 
considered as a system with n functions. The processes of changing the microroughness of the path (or changing the soil 


density) are denoted as functions q,(t) ...q¢,(t) [9, 10]. 


If the statistical description of external actions admits a stochastic orthogonal representation, then the following 


equation is valid: 
a(t)=m[aO]+ | O,(0)-4,(,0)do, 


where O;(@) is a system of random functions @; @; (t,@) 1s a system of nonrandom functions t and parameter w; m 
[g;(t)| is mathematical expectation. 

Mobile agricultural machines are systems with delayed actions [2, 6, 9], which are characterized by many 
parameters of uneven field microrelief or soil hardness [2, 8, 9]. These are one-type parameters, and they can be shifted 
in time fg. In addition, the design of the machine (front and rear axle wheels) affects the relationship between the 
actions. Given these features and the assumption that the mathematical expectation of external actions is zero, the 


matrix of spectral densities of external actions is described as follows [11, 12]: 


S,(@)  S,(@)e" C135, (J) 
S,(@, jo) =|], (@)e’" S,(@) CS, (jo)}), (1) 
C35 yp (jo) CxS, (jo) Sy (o) 


—Joto. 


= As — pJoto. = : 
5Cy, =Ley, =e" °3¢,, =1 — for rear drive 


— jot " 


where ¢,,; =t3C,,; =@ /°"3C;, =|;c,, =e/°" — for front drive axle; c,, =e 


axle. 
To simplify the calculation of statistical characteristics, we reduce the matrix to a diagonal form to obtain a 
matrix of uncorrelated external actions. We introduce the characteristic matrix and reduce the matrix (1) to the form: 
S,(@)-A(jo)  S,(@)er C,35iy, (J) 
S, (0, j@)-A( jo)-E=|| S,(@)e”  =S,(@)-A(jo) — &;5S,, (jo) |, (2) 
C3 Syn (J) CSy,(Jo) — S,(@)-A( jo) 


where 4 is a characteristic function; E is an identity matrix. 


The characteristic equation can be represented as: 





ah ( jo)+br (jo)+crA(jo)+d =0. (3) 
The characteristic functions 1, (j@),’,(j@),A;(j@) were defined using Cardano formulas [13]. 
~ 
i= The matrix of transition to a new basis constructed from the coordinates of the vectors i. ( j@) (at i=1, 2, 3) 
N 
S 
= defined from the equation (3), has the form: 
= 
§ Ni Ma Ms 
N . 
© T (jo) =WHi2 22 Nsal] - (4) 
™.v 
6. Ns 2333 
ae In this case, the inverse matrix has the form: 
47,1 J® 
(jo) = 7 ap 5) 
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where 7'( jo) is an adjoint matrix; |7'( jo)| is a matrix determinant. 


The desired matrix was determined from the expression: 
S, (jo) =T™ (jo)-S, (@, jo)-T (jo) = diag||s, (j@)-S, (j@)-S, (jo) 


where S,, S>, S3 are spectral densities of uncorrelated effects of complex-variable functions. 


(6) 








Finding the diagonal matrix (6) is simplified through solving a planar problem. In this case, external actions can 


be represented as: 


























. S,(@)  S,(@)e" 
S (@,jo)=|) ** ? ’ 7 
q (@ j@) S, (o) piers S, (o) ( ) 
Then the characteristic matrix is determined from the expression: 
. S,(@)-A(@) S,(@)e"" 
S —i(@)-E =|" | ‘ . 8 
q (o, jo) () S, (@)e” S, (@)-A(o) ( ) 
In turn, the characteristic equation can be represented as: 
A* (@)-2S, (@)-2=0. (9) 


From which, the roots of the equation are: A, (@) = 0,1, (@) = 2S, (@). 
At the root i, (o) =(), the system of equations for determining the coordinates of the eigenvectors has the 
following form: 
S,(@)n,, (j@)+5S, (@)e"" +n, (jo) =0. (10) 
Given n,, =e “°”, we find n,, (jo) =1. 


At root i, (o) = 28, (o) , we get the system of equations: 


fattest te a 


Si, (o)"” “No (jo) +); ()-N>, (jo) =0 
Then the matrix of transition to a new basis takes the following form: 


—jot (ot 
e JLo e! 0 


T, ( jo) = _] riety ‘ (12) 














Through transforming the expression (8), we find spectral densities of uncorrelated effects: 
S, (jo) =S,(@)-d+e7") =0. 
S, (jo) =S,(@)-(+e/") =0. (13) 
If we assume that the profile of the agricultural background in the longitudinal and cross section is determined 
by uncorrelated random processes [8, 9], then we can imagine the third element of the diagonal matrix as S,, (w) . In this 
case, the matrix (13) takes the form: 
S, (ja) = diag||S,(o)-(+e%");  S,(o)-(+e"); S, (o)| (14) 
Spectral densities of vibrations were determined through solving the matrix (1), together with the matrix (14) 


and with account of formulas describing uncorrelated actions. The result is the expression: 


S.(@)=S,(@)| ©}, (@)+ 7, (@)]=0. (15) 
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Fig.l Spectral velocities of vibrations in the solution using: 
curve 1— matrix (1); curve 2 — matrix (6); curve 3—matrix (15). 


Research Results. Fig. 1 shows graphs of spectral densities of vertical velocities of the front part of the 
thresher of “Vector” line harvester. Obviously, the dynamic parameters of the machine affect the vibration spectrum 
(Fig. 1). The narrowband spectrum 1s typical for conditions of close coincidence of the frequencies of natural vibrations 
of the combine axles with the reaper (Fig. 1). The width of the spectrum is also affected by the delay in the actions of 
the wheels of the combine axles. 

The results of calculating the intensity of vibrations using matrices (1), (6), (15) are diagrammatized in Fig. 2, 


which confirms the dependence of functional relationships in the machine on the nature of external actions. 


z/g 
“ETT TTT TTP Eo 
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Fig. 2. Diagrams of vibration intensity when the machine moves along microprofile with different nature of irregularities: 
1 and 2 — along the furrow; 3 — across the furrow 


Solving the systems, which are uncorrelated external actions on a moving machine with weak functional 


relationships, gives close results regardless of the methods used. 


Identical results, when using various calculation methods, are obtained only for machines with weak functional 
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relationships when moving along the field microprofile that generates an uncorrelated external action. 


The dynamic characteristics of the combine harvesters were evaluated according to the following formula: 
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W..=— | S,.(@)-do, (16) 


where w,, @> are spectrum frequency boundaries. 


The vibration intensity for various conditions determined using the formula (16), is presented in Fig. 3. 
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Fig. 3. Engineering assessment results of the harvester dynamics 


Fig. 3 shows the dependence of overloads 7, on the combine harvester velocity. Obviously, the different 


spectral composition of the irregularities of the agricultural background microprofile affects the intensity of vibrations 
of the machine nodes. There are resonant speeds in the range from 10 to 15 km/h. The graphical dependences obtained 
(Fig. 2, 3) provide estimating the load spread on the units. 
The shaded area illustrates their scattering and enables to estimate the coefficients of variation of V(z). The table 
shows parameters of the distribution of loads at medium speeds of the machines [2]. 
Table 1 


Dynamic load parameters of mobile machines of agroindustrial complex 


Fenster tice 0.840 0.332 
————e 2.781 0.835 
Tractor K-700A, 3.352 0.474 


trailer 2 PTS-4 8.323 0.510 
aus Oats casita iets [247 0.364 
es 5.527 0.382 


ee ; 2 operational (tillage) 0.752 0.346 





Results and Discussion. The analysis of the calculation results of the spectral densities of vibrations according 
to formulas (1), (6), (15) is diagrammatized in Fig. 2, 3. From the graphical dependences, it can be seen that the 
spectrum of the effects of the loads and their dynamics affect the calculation error. 

The calculation of dynamic multidimensional systems and the assessment of machine vibrations are automated 


and performed using a specialized software product. 


' Groshev LM, Partko SA, Lukonin AYu. Calculation of random vibrations of the housing of Don line combine harvester. Certificate of software 
registration no. 2012614015 Russian Federation. Applicant and proprietor, FSBEI of Higher Education “Don State Technical University”. No. 
2012611617; March 7, 2012; publ. 04.28.2012. 
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The program is used to select vibration parameters of the contours and units of mobile agricultural machines, 
to calculate the smoothness of combine harvesters, to design machine parts [14], and to assess the working conditions of 
the machine operator. 

Before using a software product, you should select a dynamic model first, then to set the design parameters of 


the machine and its units, and to determine the mass-geometric and elastic-dissipative characteristics of the system. 


Conclusions. The expressions (1), (6), (15) are acceptable for the design calculations of spectral densities of 
impacts on machine components. The use of the formula (6) is unacceptable when designing machines with a mode of 
movement along irregularities of an agricultural background with a low-frequency spectral composition and for 
machines with a low correlation of vibrations of the undercarriage and the steering. 
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Introduction. It should be noted that the study on the problem of the effect of the mass of parts on the vibration-abrasive 
processing is insufficient. In the works of A.P. Babichev and M.A. Tamarkin, the fact of such an effect 1s mentioned, 
but the degree and mechanism of the effect are not disclosed. In the metal removal formulas, only the number of 


interactions leading to microcutting is taken into account. The present work objective is to determine the effect of the 
mass of parts on the metal removal rate under vibroabrasive machining. 

Materials and Methods. An empirical, 1.e., experimental, approach is used. Parts from D16 and 30KhGSA materials which are 
widely used in the aviation industry were selected as samples. To change the mass, holes were drilled in the blanks; lead was 
poured into some samples, and plugs made of the same material as the blanks themselves were clogged into the others. Thus, 
experiments were carried out with solid, hollow, and weighted with lead samples. The working abrasive medium was scrap of 
grinding wheels of 40 x 80 mm, 25 grain size, and of trihedron prisms of 15 x 15 mm, 16 grain size. The experiments made it 
possible to clearly demonstrate the effect of grain size on the removal rate of the workpiece. 

Results. The parameters of the effect of the mass of parts on the removal rate under vibroabrasive processing are 
determined. The results obtained show the removal per unit area. The data are approximated by the least squares 
method with a linear function. A version of its distribution 1s selected using the Fisher statistical criterion. 

Discussion and Conclusion. It is shown how the workpiece mass determines the specific removal rate under the 
vibroabrasive machining. In the future, the database which is used to determine the effect of the work material 
characteristics on the process under consideration should be replenished. This will allow introducing a correction factor 
for the influence of mass in the metal removal formula, which will provide more accurate prediction of metal removal at 
the design stage of technological processes of vibration-abrasive machining. 


Keywords: vibroabrasive machining, abrasive environment, finishing and clearing treatment, mass, specific removal 
rate, approximation, Fischer criterion, influence coefficient of weight, scrap of grinding wheels, grain size. 
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Introduction. Engineering has always faced the task of improving the quality of its products. With advances in 
technology, ways of its solution are improved. This requires methods for predicting the efficiency of finishing 
processes. Vibroabrasive processing provides the required quality parameters along with high productivity, processing 
complex parts, as well as multiple-workpiece machining. In order to solve the issues under consideration, the following 





are being studied: 

- organization and development of processes and methods of affecting the working abrasive medium and the 
object to be processed; 

- development of new media and processing techniques; 

- reduction of energy costs; 

- improving the quality of processing [1]. 

During vibroabrasive treatment, metal and its oxides are removed from the surface due to mutual collisions of 
medium particles and workpieces. This process is provided through vibration of the working chamber, in which the 
workpieces and the medium are located. The camera is mounted on spring supports, so it can vibrate in various 
directions. Oscillations are transferred from an inertial (or other type) vibrator with a frequency of up to 50-100 Hz and 
an amplitude of 0.5—5.0 mm or more [2]. The number of interactions on a unit of surface of the workpiece per unit time 
is of random nature [3-— 10]. 

The study objective is to determine the effect of the mass of parts on the specific removal of metal under the 
vibroabrasive treatment. 
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Materials and Methods. The processing was carried out on a universal vibro-tumbling machine with four 
working chambers with a volume of 10 liters. The analytical balance AD 200 was used for mass measurements. 

The working abrasive medium was scrap of grinding wheels of 4080, grain size 25 [11] (Fig. 1), as well as 
trihedral prisms (TP) of 15x15 mm, grain size 16 (Fig. 2). 


rie 
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Fig. 2. Trihedral prisms, grain size 16 


After a rough turning operation, the samples were processed in an environment of trihedral prisms for 10 
minutes to remove burrs and smooth out roughness (Fig. 3). 





Fig. 3. Processed parts from aluminum D16 and steel 30GSA 


Then, processing was carried out in two stages of thirty minutes in the scrap of abrasive wheels. The working 
chamber vibrated with a frequency of 34.7 Hz and fluctuated with amplitude of 2.5 mm under continuous supply of 
process liquid (soda ash solution, 0.2%). The solution removed wear products (particles of metal and abrasive) from the 
surface of parts and the working medium. Then, processing was carried out in a TP medium (also in two stages of thirty 
minutes). 
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Research Results. As a result of the experiments, the desired values of mass m, g, and specific removal rate 
were obtained. Deviations of d from the theoretical model were estimated. The calculated and tabular values of the 
Fisher criterion were considered in comparison (Tables 1-4). 


Table 1 
Resulting data on mass and specific removal of samples from steel 30HGSA in TP environment 


p Discrepancy 8OSTEID | 
Standard deviasoon | _-401462F-06 —| 
Confidence mtonal 95% | _909679E-07 |S 
Aversa 
Right class boundary | —43i22F-O6 
/ Tet class boundary | —-2490856-06 |S 
- Feateulaed «| _——=s2.24885006—— | 
-Frabulated [a t0snossig PS 
/Angularcoeticionta | _293961F-08 |S 
/Freememberb (| —«taT9S6E-06 | 





Table 2 
Resulting data on mass and specific removal of samples from steel 30HGSA 
in wheel-scrap environment 


p Discrepancy 29RD 
/ Standard deviason | _-1s2708-06|SSS 
Confidence imenal 95% | _130868-06 | 
verge | S206 
Bish class boundary | _-6s7a7ar-o6 |S 
164 Tet class boundary | —=S.83954E-06 |S 
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Weight, m, g Removal rate, g/mm? d (deviation) 
F calculated 8.4308727 To 


F tabulated 3.105806516 re 
Angular coefficient a 4.94112E-08 re 
L97T9E-06 | 





Table 3 
Resulting data of mass and specific removal of samples from aluminum D16 in TP environment 


[Discrepancy Sd atTIEID 
Standard deviason | _1368998-06 | 
Confidence mtenal 95% | _140079F-06 | 
Avenged SS nae SSCS 
Right class boundary | —=3.5296-06 | 
Tet class boundary | 72M31E-07 | 
 Feateulaed is 93008 SSCS 
-Frabulated | szaonsssaz 
-Angularcoeticienta | _S738SF-08 | 
pFreememberb (| _—4230638-07 «YS 





Table 4 
Resulting data of mass and specific removal of samples from aluminum D16 
in wheel-scrap environment 


p Discrepancy as 
/ Standard deviason [2 798288-06 |S 
[Confidence imtenal 95% | _287576E-06 | 
averse | 9T6ASE-OS 
Rishi class boundary [‘7852198-06 |S 
ee 
ee 
ee 
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F calculated 51.65197229 





F tabulated 3.249835542 
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Weight, m, g Removal rate, g/mm* d (deviation) 


Angular coefficient a 1.37487E-07 a 
a 


We approximate the tabular data by a linear dependence using the least-squares method. We take the 





approximating function in the form: y = ax + b. Then the discrepancy (sum of squared deviations) has the form: 
S(a, b) = 3™1(y; — ax; — b)?. In the least-squares method, discrepancy should be minimal. At the minimum point of 
the multivariable function, the partial derivatives of this function with respect to independent parameters are equal to 
zero; therefore, the minimum conditions are: 


as . 
iia 0% — ax; — b)x; = 0, 
t=1 
n 
OS 
aa 22.01 —ax; —b) =0. 
i= 


After transformations, we obtain the following system of two algebraic equations with two unknowns: 
{" Dit, PD Liar ke = Dixy (1) 
adie +bn =D. 

Denote the mass of parts by x, y is the specific removal of workpieces. We approximate the given tabular 
function by a linear dependence. To determine the best parameters a and b by the least-squares method, we solve the 
system (1). We solve the system using the matrix method in Microsoft Excel and obtain the values of a and 5 (see 
Tables 1- 4). 

To check the adequacy of the results, we use the Fisher criterion and tabulate them (see Tab. 1-4). The 
calculated value of the Fisher criterion has the form: 

(Vi pacuer. ~ Yepenuee pacuer.)~ n-t-1 
Moacuer. = . a * X(Vi - Vipacuer.)” 
where ¢ is the number of factors x affecting y; n is the number of observations. 

Through comparing the calculated and tabulated values of the Fisher coefficient (see Tab. 1-4), we see that the 
calculated F significantly exceeds the tabulated F’. Thus, we can conclude that the constructed dependence corresponds 
to the initial data with 95% confidence. 

Substitute the obtained values of a and 5. 

For 30KhGSA samples in the TP environment, y = 2.93961E — 08 x x + 1.47936E — 06, in the wheel-scrap 
environment, y = 4.94112E — 08 xX x + 1.97429E — 06. 

For Dl6samples in the TP environment, y = 5.73851E — 08 xX x + 4.23063E — 07, in the wheel-scrap 
environment, y = 1.37487E — 07 xX x + 9.04509E — 07. 

Graphically, the results are presented in Fig. 4—7. 
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Fig. 4. Graph of dependence on specific removal mass of material of 30HGSA samples processed in TP 
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Fig. 5. Graph of dependence on specific removal mass of material of 30HGSA samples processed in wheel scrap 
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Fig. 6. Graph of dependence on specific removal mass of material of D 16 samples processed in TP 
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Fig. 7. Graph of dependence on specific removal mass of material of D 16 samples processed in wheel scrap 
The graphs show that the specific removal rate varies significantly depending on the environment, as well as 
on the mass and material of the workpieces. 
Conclusion. In analyzing Fig. 4-7, we can conclude that under machining workpieces of a larger mass, 
specific removal rate increases. This is due to changes in the momentum. 
Since the momentum is equal to the product of the mass of the body and its speed, then with an increase in the 
mass of the workpiece, the momentum of the interaction of particles and the workpiece surface increases. With an 


increase in the mass of workpieces twofold, the specific removal rate increases by 1.5—2.2 times. When comparing the 
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values of the angular coefficient a and the free term b between treatments, we can assume that they are affected by the 
parameters of the graininess of the media, as well as the characteristics of the processed material. The results obtained 
allow us to verify the theoretical dependences through introducing a coefficient to determine the impact of mass ratios. 
This will provide more accurate prediction of removal at the design stage of the vibration-abrasive processing. The 
results obtained replenish the database, which is used to determine the effect of the work material characteristics and the 


environment on the process under consideration. 


References 

1. Babichev AP, Babichev IA. Osnovy vibratsionnoi tekhnologii [Fundamentals of vibration technology]. 2nd 
ed. Rostov-on-Don: DSTU Publ. House; 2008. 3 p. (in Russ.) 

2. Babichev AP, Motrenko PD, Ivanov VV, et al. Vibratstonnaya mekhanokhimiya v protsessakh otdelochno- 
uprochnyayushchei obrabotki 1 pokryti1 detalei [Vibration mechanochemistry in the finishing and hardening processing 
and coating of parts]. Rostov-on-Don: DSTU Publ. House; 2012.5 p. (in Russ.) 

3. Martynov AN. Osnovy metoda obrabotki detalei svobodnym abrazivom, uplotnennym inertsionnymi silami 
[Fundamentals of the method of machining parts with a free abrasive compacted by inertial forces]. Saratov: Izd-vo 
Saratov. un-ta; 1981. 212 p. (In Russ.) 

4. Babichev AP, Ryseva TN, Samadurov VA, et al. Naladka 1 ehkspluatatsiya stankov dlya vibratsionnoi 
obrabotki [Adjustment and operation of machines for vibration processing]. Moscow: Mashinostroenie; 1988. P. 12. (in 
Russ.) 

5. Tamarkin MA, Krashenitsa SB. Vliyanie mikrorel'efa chastits rabochei sredy na_ proizvoditel'nost' 
vibroabrazivnoi obrabotki [Influence of the microrelief of working medium particles on vibroabrasive machining 
efficiency]. In: Progressive Hardening Technology. Rostov-on-Don: DSTU Publ. House; 1981. P. 12—15. (In Russ.) 

6. Tamarkin MA. The optimization of technological processes of details processing by free abrasives. Key 
Engineering Materials. 2005;291/292:3 19-322. 

7. Tamarkin MA, Tishchenko EE, Rozhnenko OA. Metal removal in the abrasive machining of complex 
surfaces. Russian Engineering Research. 2013;33(5):302—305. 

8. Clark J, Massarsky ML, Davidson DA. Edge and Surface Conditioning for Improved Part Performance and 
Service Life. Products Finishing. November, 2016. P. 2. 

9. Massarsky ML, Davidson DA. Turbo-Abrasive Machining and Finishing. Manufacturing Engineering. June, 
2014. P. 22. 

10. Massarsky ML, Davidson DA. Free Abrasives Flow for Automated Finishing. Manufacturing Engineering. 
May, 2013. P. 43. 

11. Kashchuk VA, Vereshchagin AV. Spravochnik shlifovshchika [Handbook of the grinder]. Moscow: 
Mashinostroenie; 1988. P. 19. (In Russ.) 


Submitted 09.04.2020 
Scheduled in the issue 29.05.2020 


About the authors: 
Butenko, Viktor I., professor of the Engineering Technology Department, Don State Technical University (1, 
Gagarin sq., Rostov-on-Don, 344000, RF), Dr.Sci. (Eng.), professor, ScopusID 35972120700, ORCID: 


https://orcid.org/0000-0001-9260-1030, butenkowiktor@yandex.ru. 


Butenko Vv. L, et al. Effect of mass of parts on removal rate under vibroabrasive machining 





Stel'makh, Andrei V., graduate student of the Engineering Technology Department, Don State Technical 
University (1, Gagarin sq., Rostov-on-Don, 344000, RF), ORCID: https://orcid.org/0000-0002-8630-8748, 


andreystelmakh9 1 @mail.ru. 


Claimed contributorship 
V. I. Butenko: academic advising; analysis of the research results; the text revision; correction of the 
conclusions. A. V. Stel’makh: basic concept formulation, research objectives and tasks setting; computational analysis; 


text preparation; formulation of conclusions. 


All authors have read and approved the final manuscript. 


Machine building and machine science 


—_" 
ON 
\O 


Vestnik of Don State Technical University. 2020. Vol. 20, no. 2, pp. 170-177. ISSN 1992-5980 eISSN 1992-6006 











MACHINE BUILDING AND MACHINE SCIENCE 


UDC 621.791 


Analytical determination of wear resistance criteria of a stamping 
tool for various loading conditions 


S. D. Kolotienko, A.V. Zhuravlev, E.V. Roshchina 


Don State Technical University (Rostov-on-Don, Russian Federation) 


Introduction. The paper is devoted to the analytical determination of general criteria for predicting the wear resistance 
of heavily loaded friction couples, in particular, stamping tools, operating under various loading conditions. The 
proposed system of criteria is based on physical dependencies that link the basic wear index, i.e. the number of loading 
cycles N; with crack resistance parameters n and C. The work objectives were the development and analytical 
foundation of the calculation technique that provides for predicting the wear resistance of the stamping tool for various 
loading conditions, as well as the predetermination of experimental verification of the efficiency of the developed 
methodology. 

Materials and Methods. Mathematical models that link the key criterion of wear-resistance, the number of loading 
cycles N;- with crack resistance parameters n and C, are proposed. An analytical verification of the proposed models is 
carried out. 

Results. Mathematical models are developed for predicting the wear resistance of a stamping tool operating under 
various loading conditions. In particular, for the following cases: under sliding and rolling friction of the tool on a 
plastically deformable metal, for conditions of thermo-mechanical contact fatigue, for high-cycle brittle contact 
damaging, for pulsating contact on a plastically deformable metal, as well as for combined types of frictional contact. 
Discussion and Conclusions. The results obtained can be used in the design and structural optimization of the stamping 
tool operating under various loading conditions, as well as in predicting its life cycle. 


Keywords: stamping tool, wear resistance criteria, loading conditions, contact damaging, crack resistance, sliding, 
rolling, surfacing materials. 
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Introduction. Experimental and theoretical studies on the wear processes of surfacing materials show that no 
material properties separately can determine the wear resistance under the frictional contact extreme conditions in a 
unique manner. Only a certain set of physical and mechanical properties with account of the frictional interaction 
features should be considered when determining the wear resistance of surfacing materials. 

Wear resistance criteria should meet the conditions of specific tribocouplings. To determine the criteria, a 
single approach should be followed in each case. Obviously, for all the phenomena and types of contact damage, the 
principle condition should be taken into account: that material is considered more wear-resistant whose stage of fatigue 
crack growth is longer. Thus, the wear resistance criterion should take into account the duration of the process of 
destruction of the surface layer and separation of the first wear particle. In this case, the determination of the wear 
resistance criterion is reduced to evaluation of the cycle index N; at which the surface layer destruction occurs. 

In this case, the assignment of wear resistance criteria may consist not in calculating the contact endurance of 
materials, but in comparing and selecting surfacing materials for specific frictional contact conditions that are 
characteristic of heavy loaded friction couples. Thus, the criteria should provide adequacy of the wear resistance series 
for specific types of friction contacts. 

In accordance with a unified approach, the following key points should be taken into account when 
determining the criteria. One type of crack should be considered in all types of frictional contact. Generically, the crack 
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resistance parameters of materials that determine the kinetics of the development of fatigue microcracks should be taken 
into account. When comparing surfacing materials with each other according to the criteria, it is required to consider the 
170 possibility of developing a competing type of wear — low-cycle viscous damage as a result of the petal peeling. 
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Materials and Methods. For reversible sliding conditions on a plastically deformable metal, the author first 
proposed a method for analyzing wear resistance considering the basic properties of the material, which determine the 
surface layer destruction under friction [1]. 

It appears possible to determine the wear resistance criterion with account to the fracture toughness, yield 
strength, friction coefficient, and initial crack sizes. The studies [2—9, 11, 13, 15, 16] show that this approach is also 
promising for other frictional contact conditions characteristic of heavy loaded friction couples including a stamping 
tool. 

V.V. Rubanov, assuming that the criterion proposed for reverse sliding in [1], and using the dependence of the 
yield strength on a number of characteristics of the structure of the material, proposed to determine the wear resistance 
of the material according to the criterion [4]: 


‘6. +11Gb +) -K.. 
K, = a) ” 
where o, is the ultimate mechanical stress in the material from the action of the matrix under which its destruction 


occurs; /; is the average defect size; K,, is fracture toughness; @ is the volume fraction of defects; G is the shear 
modulus; b is Burgers vector; fis the friction factor (from now on, the friction factor is applied on the ground of the 
Siebel law). 


Given that for a number of alloys, there is the dependence K,, on the tensile strength under bending o _ , after 


a series of simplifications, V.V. Rubanov recommends the following criterion for assessing the wear resistance of 
surfacing materials: 
G:0:0 
Eg as : aT 
f° --9)-f, 


This approach has limited possibilities due to the following. The proposed criterion is common for various 


(2) 


frictional contact schemes; and when determining the wear resistance of surfacing materials for sliding and rolling 
according to the criterion (2), differences in the contact damage mechanism are not taken into account under the 
conditions of contact interaction for sliding and rolling, both with tangential with effort or without it. 

The wear resistance criterion (2) does not take into account the crack resistance parameters of materials that 
affect directly the rate of development of fatigue microcracks, and, therefore, determine the duration of contact damage 
processes. 

Research Results. In the analytical definition of wear resistance criteria, the basic assumption should be made: 
the field of contact tensile stresses is uniform, and its value is equal to the maximum value that is determined on the 
surface. Thus, this assumption does not take into account the gradual attenuation of tensile stresses in the depth of the 
material and is a more stringent estimate. The studies [10] show that the errors arising in this case are small if the length 
of the cracks is less than the depth of the zone of compressive stresses. 

The task of determining the criteria for the wear resistance of materials wearing under the contact damage 
conditions is as follows. It is required to analytically determine the time (the number of contact interaction cycles) of 
separating the particles of abnormal wear. In this regard, the stresses that act on the crack in a particular frictional 
contact should be used as tensile stresses. 

Given the differences between low-cycle brittle contact damage and multi-cycle damage, the problem of 
determining the criteria for wear resistance is divided into two types. In the first case, when determining the contact 
fracture process time for low-cycle brittle contact damage, it is required to consider the possibility of reaching a critical 
crack. In the second case for low-cycle brittle contact damage, it was found that cracks do not reach critical values 
(/.>>1;). Here, it is possible to use the maximum permissible size of spalling as the final crack length (it is established 
either from the requirements for wear tolerance of tribocoupling or from the requirements for limiting the size of 
contact damage on the edge of the working surface). 

Wear resistance criterion for sliding on plastically deformable metal 

The contact fracture process time is determined through integrating the Paris equation [17]. As a value of 
tensile stresses for this case of frictional contact, we should take [10]: 
_40f 


pmax 
TU 


(3) 


where p is the average pressure in contact. 
In this case, the expression for the stress intensity factor takes the form: 
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The fracture process time, and, consequently, the wear resistance criterion N;, is defined as 








kdl wk dl 
Nae a) (5) 
1, CK" C(4pf)i, 
At n# 2, the expression (5) solution takes the form 
2n"* l l 
n—2 |° (6) 


(n—2)C(4pf) 12 72. 


From the basic principles of fracture mechanics, it is known that the critical crack length is the value at which 
the crack starts a supercritical self-motion. This moment determines uniquely the fracture toughness of the material Kjc. 
In the case of contact loading by a uniform stress field under sliding, the value of the critical crack length can be 
determined from the expression [10]: 








TU f" 
= lop p (7) 
After the corresponding replacement and transformations, we obtain 
Noe an | tf _4of - . (8) 
© (n= 2)CA4PPY" | A (ePKic 


This expression can be used as a criterion of wear resistance when sliding against a plastically deformable 
metal under the conditions of high contact pressures, when low-cycle brittle contact damage is realized. In this case, as 
follows from the expression (7), the wear resistance of surfacing materials is determined by the contact pressure, 
friction coefficient, fracture toughness, size of structural defects and crack resistance parameters n and C. 

Wear resistance criterion for thermomechanical contact fatigue 

The effect of heating-cooling cycles, which is characteristic primarily for a pulsating contact with a heated 
plastically deformable metal, on the friction surface causes the origination of periodic thermomechanical stresses in the 
surface layer. The expression of V. I. Dukhovchenko [12] for determining this type of stress 1s known as follows: 

_ aEAT (9) 


l-yv 





where « is the thermal coefficient of linear expansion; E is the elasticity modulus; AT’ is the temperature difference of 
the heating-cooling cycle; v is the Poisson's ratio. 
In this case, with regard to the expression (9), the solution to the integral for the Paris equation at n# 2 is 


determined from the formula: 
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If the limiting value 5,,,,., which is assigned as the tolerance on the limiting size of the centre of spalling of the 
working surface edges, is taken as /., we obtain 
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This expression can be used as a criterion for the wear resistance of surfacing materials under 
thermomechanical fatigue. Under these conditions, wear resistance of materials also depends on the crack resistance 
parameters n and C, is determined by the tolerance for unit spalling, the elasticity modulus, the Poisson's ratio, and is 
inversely proportional to the coefficient of thermal expansion, the magnitude of the temperature difference and the size 
of the initial defects. 
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Wear resistance criteria for multi-cycle brittle contact damage 

We consider /,= /; and take into account that when rolling with a tangential force and an elastic pulsating 
contact, the angle of crack development does not depend on the material properties and is close to 45°. Then, the 
following expression is obtained to determine the fracture process time under rolling with a tangential force: 























ae: 7.35 | 1 1 13) 
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For rolling without tangential force 
2 |4627} 1 1 
N= | 42 ma ed |? (14) 
(n—2)C 0 , max 12 ] 2. 
Upon the pulsating contact 
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Since, under conditions of multicycle brittle contact damage, the cracks do not reach a critical value during 
development, /.>>/;, the expressions (13) - (15) can be transformed as follows: 


— for tangential force rolling 








21.30" 
N7= ee (16) 
(n—2)CO+IISY Grant; ” 
— for rolling without tangential force 
2-4.62” 
$s (17) 
(n 2)COe max! : 
— for pulsating contact 
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The obtained criteria for wear resistance indicate that in this case, too, the critical impact is exerted by the 
parameters of crack resistance of materials n and C, the size of the defects; and, in addition, in the case of rolling with a 
tangential force, the friction coefficient should be taken into account. 

Wear resistance criterion for rolling against plastically deformable metal 

A feature of this frictional contact type is that the counterbody acquires intense plastic deformation. In this 
case, stresses in the contact zone cannot be estimated through the formulas for elastic contact, and the wear resistance 
criteria (16) - (18) are not applicable. 

To solve this problem, we use the analogy for rolling against a plastically deformable metal and sliding under 
the same conditions. Then, the crack will be periodically affected by a plastically deformable rolling body. Considering 
that rolling under these conditions is accompanied by significant spalling, we use the expression (6) for the sliding 
contact. As the final crack length, the limit tolerance for unit spalling of the working surface edge should be used. 

In this case, the number of cycles for the destruction of the surface layer can be determined through the 











expression: 
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After transformations, the expression (19) can be rewritten as follows: 
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Thus, the wear resistance of surfacing materials when rolling against a plastically deformable metal is 
determined through the crack resistance parameters n and C, contact pressure, friction coefficient, limit tolerance on the 
size of the spalling and the sizes of the initial defects. 

Wear resistance criterion for pulsating contact on plastically deformable metal 

The characteristic features of this type of friction contact are extremely high contact pressures (up to 4000 
MPa), which is associated with the dynamic nature of the application of the contact load. In this case, making use of the 
expression (18) for an elastic pulsating contact is also impossible. 

For the maximum tensile stresses acting on the crack, we take the value that, at a first approximation, can be 


determined using the condition o_.. ~p from the formula: 


Z Max 


1-2v 
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where p is average contact pressure. 
The fracture process time is determined through the integration of the Paris equation for the given conditions of 


contact loading: 
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Let us present the value of the critical crack length /, through K/ : 
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3 
and we get after the transformations: 
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The given expression is a criterion for the wear resistance of surfacing materials for pulsating contact against 
plastically deformable metal. Under these conditions, wear resistance is determined by fracture toughness, crack 
resistance parameters, Poisson's ratio, contact pressure and the size of the initial defects. 

Wear resistance criteria for integrated friction contact modification 

Under the conditions of sliding friction against a heated plastically deformable metal, the tool working surface 
is subjected to alternating stresses from friction forces and thermomechanical stresses. In this case, according to the 
hypothesis of a linear summation of stresses, the stress intensity factor considering the expressions (3) and (4) can be 
represented as: 

K =o Vnd =( SPE, SEAT) aj (25) 
TT l-v 
Then, through integrating the Paris equation, we can obtain an expression for determining the number of cycles 


necessary for breaking: 
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The critical crack length can be determined from the expression (25) given that the supercritical fracture point 


occurs at K = Kj. Then 
K/ 
1 =———_—!__.. (27) 
“(“ Pf, ear | 
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Substituting the value /, into the equation (26), we can obtain an expression for the wear resistance criterion in 


the case of a combined effect on the sliding friction surface and thermomechanical stresses: 
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A combined effect on the friction surface is possible under an elastic pulsating contact with the application of 


thermomechanical stresses. In this case, the equation (18) should be written as: 
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Given that in this load case /,.>>/;, we obtain the expression for the wear resistance criterion for elastic 


pulsating contact with thermomechanical fatigue 
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Discussion and Conclusions. Analytical dependences are obtained for calculating the wear resistance index 
N- for various working conditions of a stamping tool based on general criteria describing the crack resistance 
parameters of materials n and C. The results can be used in the design and optimization of a stamping tool operating 
under various loading conditions, as well as under predicting its life cycle. 
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Algorithm of software package of intellectual decision support when 
designing cyber security system at the enterprise 


E. A. Vitenburg, A. V. Nikishova 
Volgograd State University (Volgograd, Russian Federation) 





Introduction. To increase the decision-making efficiency at the enterprise, it is advisable to use a special software pack- 
age of intellectual support. Such a product is necessary when designing an information security system and increasing 
its invulnerability during modernization or configuration changes. Research objectives are as follows: to develop an 
algorithm and a mathematical model of the software package for intellectual decision support. 

Materials and Methods. The decision support method under designing an information security system is based on the 
use of a neural network (multilayer perceptron). For an objective assessment of the initial security of an information 
system (IS), a mathematical model for the analysis of security events is developed. 

Results. The statistics of malicious attacks on the IS of enterprises is analyzed. The need for timely and accurate mod- 
ernization of the information protection system is determined. Important characteristics of the designing an information 
security system are the speed at which the result is obtained and the reduction in the residual risk of IS. In this regard, 
the use of artificial intelligence systems in the process of determining the best set of protection subsystems is important. 
The threats to cyber security (CS) are classified. The main classes of security events are defined. A mathematical model 
of the neural network is developed; the input parameters of its operation are indicated. The current enterprise IS gener- 
ates numerous events which necessitates the automatic collection and analysis of data from subsystems for registering 
IS objects. The process of analyzing security events is considered in detail since the adequacy of the generated design 
decisions depends on the correctness of the data obtained in this way. The algorithm of the software package is formed. 
Discussion and Conclusions. The results can be used in the design of the information security system at the enterprise. 
In addition, CS administrators can use the developed software package to adjust the configuration settings of 
information security tools. The proposed solution will minimize the destructive influence of the developer of the 


security system which may and happen to be subjective. 
Keywords: cyber security, information system, neural network engineering, multilayer perceptron, algorithm. 
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Introduction. With the development of industry in Russia, the number of enterprises classified as objects of 
critical information infrastructure (OCII) is growing. Statistics on the distribution of attacks from the Positive Technol- 
ogies' vendor shows an increase in the number of successfully implemented malicious actions in this area. In 2019, 125 
attacks on industrial information systems (IS) were recorded. This is more than three times (or by 212%) higher than the 
same indicator in 2018 (40 attacks). The diagram of the distribution of the number of attacks by quarters of the year (Q) 


is shown in Fig. 1. 
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Fig. 1. Number of attacks on IS of enterprises in 2018 and 2019 


Industrial IS are attacked mainly using malicious software (90% of attacks). This area is supervised by 
the Federal Security Service and the Federal Service for Technical and Expert Control. The need to ensure 
cyber security (CS) of industrial information systems is confirmed by statistics of the head research centers [1]. 
According to the regulatory legal acts of the Russian Federation’, the owners of the OCII use a set of organiza- 
tional and technical measures to ensure the safe operation of the information infrastructure. At the same time, 
the legislation provides for a periodic audit of the information protection system (IPS) performance, an assess- 
ment of its effectiveness. This results in prompt updating of the configuration settings of available tools or the 
retrofitting of the system with information security tools. A special role is played by the efficiency and accura- 
cy of the decisions made, and the value of residual risk should not exceed the established indicators [2]. In this 
regard, it is proposed to automate the decision support process under designing the IPS at the enterprise. 

Materials and Methods. To support decision-making in the design of IPS, a method based on a neural 
network (multilayer perceptron) was used [3]. Input data for the operation of a neural network were information 
security threats and security events. In addition, for an objective assessment of the initial information system 


security, a mathematical model of the analysis of security events was developed. Within this model, the calcu- 


'Aktual'nye kiberugrozy: IV kvartal 2019 goda [Current cyber threats: Q4, 2019]. Positive Technologies. URL: https://www.ptsecurity.com/ru- 
ru/research/analytics/cybersecurity-threatscape-2019-q4 (accessed 24.02.2020). 

*O bezopasnosti kriticheskoi informatsionnoi infrastruktury RF: feder. zakon ot 26.07.2017, Ne 187-FZ [ On security of critical information infrastruc- 
ture of the Russian Federation: Federal Law of July 26, 2017, no. 187-FZ]. RF State Duma, Federation Council. FSTEC of Russia. URL: 
https://fstec.ru/tekhnicheskaya-zashchita-informatsii/obespechenie-bezopasnosti-kriticheskoj-informatsionnoj-infrastruktury/285-zakony/1610- 
federalnyj-zakon-ot-26-1yulya-2017-g-n-187-fz (accessed 18.05.20). 
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lated measures of IS similarity are compared with one of the IS security levels. Weighted Manhattan Distance 
is used as a similarity metric. 

Study Results. In the framework of this study, an approach to the design of a protection system based 
on the factors of the cyber security subsystem importance included in the IPS [4] is proposed. With these data, 
it is suggested to make a list of security tools for the most critical subsystems. This approach provides strength- 
ening the IPS through neutralizing critical threats. Moreover, criticality should be specified through an auto- 
mated analysis of security events [5, 6]. The vector of cyber security subsystem importance is generated ac- 
cording to the expression 

V = S(Class_Thr), (1) 
where V = (V,,...,V,) 1s the vector of IPS subsystem importance; Class_Thr = (Class_Thr,,...,Class_Thr¢) is 
the vector of criticality of threat classes; S is the functional relationship defined by a neural network. 

For convenience, we divide set of CS threats Thr into classes Class_Thr, which are defined as in [6]: 

Class_Thr = { Thr_Br,Thr_L,Th_Dist, Thr_Loss, Thr_B,Thr_A}, (2) 
where Thr_Br is the class of CS threats of the “cracking” type; Thr_L is the class of CS threats of the “leak” 
type; Th_Dist the class of CS threats of the “distortion” type; Thr_Loss is the class of CS threats of the “loss” 
type; Thr_B is the class of CS threats of the “blocking” type; Thr_A is the class of CS threats of the “abuse” 
type. 

To generate the vector of IS CS threat criticality, it is proposed to form two matrices of compliance: 

— set of threats to set of classes of threats — MThr, 

— set of security events to threat classes — MEvent. 

MThr matrix: 

MThr = Thr x Class_Thr = (mth,;). (3) 


Here, mth,; is determined from the formula: 


i 
mth, = fr if the threat belongs to the threat class, (4) 
O if not. 
Here, 7 1s the number of threat classes to which the threat belongs. 
In this case, Vi, )}; mth,; = 1. 


MEvent matrix: 


MEvent = Events xX CThr = (mevj;;). (5) 

Here, mev;; is determined from the formula: 
_ ¢1,if an event occurs when a class threat is implemented, 6 
ie aoe 0 if not. (©) 


To form the vector of the cyber security subsystem importance, it is proposed to use a neural network in a 
software package — a multilayer perceptron operating according to the formula given in [7]: 


INox = Vx 
Outi; = fQiwiylnij — 94), (7) 
Injj = Outj_1) 


where [Nox is the k-th neuron of the input layer; v, is the k-th element of the input vector; Out;; is the output value of 
the j-th neuron of the i-th layer; f is the neuron activation function, which is determined by the functional dependence 
V (1); wij, is the weight of the /-th input of the j-th neuron of the i-th layer; Inj; is the value of the /-th input of the j-th 
neuron of the i-th layer; 0;; is the activation level of the j-th neuron of the i-th layer; Out,_1; is the output value of the /- 


th neuron of the (i — 1)-th layer. 
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At the stage of analyzing security events, the initial data are the security event logs created by the system and 
the application software of the enterprise IS. Set of the security Events include the following classes [8]: 
Events = {EnterEv, ManagementSubEv , AccessObjEv , PolicyChangeEv , 
Events = {EnterEv, ManagementSubEv , AccessObjEv , PolicyChangeEv , 
UsePrivilegesEv, ISProcessesEv, LevellSEv}. (8) 
Here, EnterEv are events of the “subjects’ log-on” class; ManagementSubEv are events of the “subject manage- 
ment” class; AccessObjEv are events of the “accessing objects” class; PolicyChangeEv are events of the “system 
policy change” class; UsePrivilegesEv are events of the “subject’s using exclusive privileges” class; [SProcessesEv 
are events of the “system processes running” class; LevellSEv are events of the “system level” class. 
Dangerous events are selected from the set. Matrices of compliance are used to compare the set of dangerous 


events to threat classes. The matrix developed in the framework of this study is presented in Table 1. 


Table 1 
Matrix of compliance of CS threats of enterprise IS 


Threat types 


Information technology, computer science, and management 


Security event class 


ManagementSubEv 
PolicyChangeEv 
UsePrivilegesEv 
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posed to the urgent threat classes. Compliance is determined using a neural network with account for the functional 
dependence V in accordance with the formula (1). 

Quite an important factor is the number of analyzed sets of events and their sources. It should be noted 
that the number of security events is directly proportional to the number of sources — information resources of 
the enterprise IS [8]. Given a large number of events generated by a working IS, it is valid to automatically col- 
lect and analyze data from the registration subsystems of IS objects that describe events. In this regard, it is 
worth considering in detail the analysis of security events since the adequacy of the generated design decisions 
depends on the correctness of the conclusions received at this stage. 

According to [9, 10], at a given instant 7, the current state of the enterprise IS S; € State, State = 
{Snorm, Sdang, Sanorm} can be characterized as: 

— normal (Snorm) — normal system operation in accordance with its tasks and as given in the docu- 
ments regulating the work; 

— dangerous (Sdang) — incorrect IS operation, when malfunctions associated with hacker attacks, 
crashes and failures of software and (or) hardware are recorded; 

— abnormal (Sanorm) — temporary change in the normal IS operating mode and a surge in abnormal 


activity of users, programs and network traffic. 
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Dangerous and abnormal events require a more detailed analysis. A set of such events is the input to the stage 
of matching threat classes. Such events indicate the implementation of security threats. 

Therefore, the results of monitoring and analysis of security events will be the input to the neural network (7). 

Any EventIS; security event can be described by an attribute tuple [9]: 

EventIS; = (ID, Data, Level, Source, EventT ype, EventState, SecureParams), (9) 
rae JD is the event ID; Data is the event generation time; Level is hazard level of the event; Source is the 
source of the event; EventType is the event type; EventState is event status; SecureParams is event security 
parameter vector. 

SecureParams = (h,u, risk), (10) 
where / is an indicator for the generation of the event of a certain code relative to the total number of security 
events over the period A7; u is the severity of the consequences of the event (potential damage); risk is the IS 
information security risk. 

Over the period AT, a set of IS events (Event/S) which should be evaluated to determine the level of IS 
security is generated in the information system. 
The ratio of the number of events of one type or another to the total number is determined by the indi- 


cator h: 


NEventID 
h=—_ (11) 
NEvent 


where NEvenID is the number of events of a certain code, NEvent is the total number of events over the period 
AT. 

The application software defines event codes differently. So, in the framework of this work, we have 
considered the encoding of the Windows OS Microsoft’. 

The IS information security risk 1s a function of the frequency of the event and the potential damage: 

risk =hxX wu. (12) 

The sum of private indicators determines the overall risk RiskSum: 

RiskSum = Vive” risk;. (13) 

The set that provides the classification of events and the assessment of private indicators of IS security 
is described by subsets of elements: 

PPS = {{EvType}, {EvState}, {ISState}}. (14) 
Here, EvType is a set of types of the events detected; EvState is a set of possible event states: 

EvState = {Ev", Ev’, Ev}, (15) 
where Ev™ — events of normal IS operation; Ev“ — IS CS threat events; Ev“— abnormal events characterizing 
deviations of the IS from the normal mode of operation (additional analysis is needed). 

To determine whether event Event/S; belongs to one of the three states EvState, the classification prob- 
lem is solved — the set SP = DamgeEv U NormEv 1s used, which is divided into two basic subsets. They are 
formed by an expert group on the basis of data on quick events of the normal IS operating mode and the previ- 
ously detected CS attacks and incidents that describe the signature of IS quick and threat events: 

— DamgeEv is a set of events that are known features of an attack or determine the scenario of an in- 
cident; 


— NormEv is a set of events typical for the normal mode of IS operation. 


*Opisanie sobytii sistemy bezopasnosti v Windows 7 i Windows Server 2008 R2 [Description of Security Events in Windows 7 and Windows Server 
2008 R2]. Microsoft. URL: https://support.microsoft.com/ru-ru/help/9775 19/description-of-security-events-in-windows-7-and-in-windows-server- 
2008 (accessed 05.08.2019). 
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F(Event!S;(EvState) — {SP},;) are classified as follows: 


Ev", EventlS; € NormEv, 
Event!S,(EvState) = Ev‘, EventlS; € DamageEv, (16) 
Ev*, EventlIS,; € NormEv U DamageEv. 


Here, EventlS;(EvState) are IS events with a status attribute, each of which corresponds to a set of connections 
SP; with quick events from a set of templates of the normal mode and the mode with CS breach. 

If the event is not present in the profiles of quick events or security breach events, it 1s determined to 
be abnormal. The reasons for its occurrence should be considered separately by the administrator of the IS sys- 
tem. 

ISState ={ISnorm, ISdang,ISanorm} is a set of IS statuses indicating a normal mode of IS operation. 
The IS status is determined from the formula: 


ISnorm, VEventlS; € EventIS|EvState = Ev", 
ISState = { ISdang, JEventIS, € EventIS|EvState; = Ev®, (17) 
ISanorm, ISnorm 0 ISdang, 


where /Snorm is the mode of normal IS operation; JSanorm is the mode of IS operation with features of abnor- 


mal activity; [Sdang is the mode with fixed IS security breaches (Fig. 2). 


Described by set 
Described by set of security 
of normal events ISdang events 
TSanorm 
Described by set 
of abnormal 
events 


Fig. 2. Relationships of subsets of IS states 


Based on the data obtained, a vector is formed that determines the IS status, the total risk, as well as 


the proportion of abnormal and security breaches: 


IS = USState, RISKsum, NanalEv, NdandEv), (18) 
where NanalEv is the proportion of detected abnormal events, NdandEv is the proportion of events of CS 
breach of the enterprise IS. 

To make a decision on whether to strengthen information security subsystems, the IS security level SZ 
is calculated: 

SL = {safe, stable, abnormal, crisis, dangerous}. 

To determine the level of IS security, the designer of the CS system generates a vector that defines the 
reference indicator of the IS security [S_Perf. SL is formed according to the similarity of vectors JS and 
IS_Perf. 

Five levels are accepted for evaluating the IS security within the framework of this mathematical mod- 
el; therefore, when determining whether the IS belongs to one of the levels, the k-nearest neighbor method is 


used [10]. To this end, intermediate vectors corresponding to the remaining four security levels are compiled. 
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This operation is based on the values of the /S_Perf vector corresponding to the safe level of the IS status when 
multiplied by a scalar correction factor (the value is determined by experiment): 
kx RIS = (k x [SState,k x RISKsum,k x NanalEvent,k xx NdandEvent). (19) 
Weighted Manhattan distance 1s used as a similarity metric [10]: 
p(IS, RIS) = w Y4_,|1S; — RIS)|. (20) 
Based on the presented mathematical model, a generalized algorithm for the operation of a software 


prototype has been developed (Fig. 3). 
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Fig. 3. Algorithm for software prototype of decision smartness complex when designing an information 
security system at the enterprise 
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Step 1 (blocks 1, 2). Launching a software prototype. The algorithm start. Input of data on the require- 
ments vectors /S_Perf to the IS status of each security level. Enter the period over which monitoring will be 
conducted. Loading from the DB templates with sets of normal and dangerous events. 

Step 2 (blocks 3—5). Monitoring OS event logs, receiving records on each event. Event listing. 

Step 3 (block 6). Data analysis of the collected events, formation of the event tuple — the formula (8). 
Classification of events according to the formulas (14—15) into normal, abnormal, and dangerous ones. Calcula- 
tion of the frequency of occurrence of events, potential damage and risk using the formulas (11-13). 

Step 4 (blocks 7-11). Classification of IS states according to the formula (17) based on the event dis- 
tribution data into a set of abnormal, normal, and dangerous ones (step 3). 

Step 5 (blocks 12—14). Based on the calculated data and event tuple in the IS, the formation of the cur- 
rent status vector, calculation of similarity between JS and [S_Perf vectors: the formulas (19), (20). Decision- 
making on whether IS belongs to one of the five security levels. 

Step 6 (block 15). The algorithm completion. 

Discussion and Conclusion. In the framework of this study, an approach to modeling an information 
security system is proposed. Different qualitative and quantitative complex of security tools is taken into ac- 
count depending on the actual threats to information security breaches. The presented method provides increas- 
ing the efficiency of the introduced cyber security system and reduces the likelihood of a designer mistake. The 
software package developed as part of this study has advantages that cannot be obtained under “manual” de- 
sign: 

- accounting of all data on the protected system, 

- obtaining accurate results at the earliest. 
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Introduction. Digital systems that control the maintenance of separate mechatronic process facilities (MPF) and sets of 
production machines are mainly considered. Numerous issues on maintaining the reliability of the condition and 
emerging malfunctions, as well as the multifactorial nature of the using the existing monitoring and diagnostic systems, 
are noted. In this regard, the relevance of the tasks of developing methods of processing equipment maintenance to 
make decisions under the data veracity and limitation is specified. 

Materials and Methods. To analyze the criticality of the technical condition, an assessment of the efficiency of the 
autonomous control of the device state 1s formed. The method of the neuro-fuzzy system is used to determine the 
aggregate criterion of criticality. It 1s proposed to apply this approach to develop recommendations on equipping a 
production facility with the necessary means of maintaining overall performance and reliability. 

Results. The solution provides predicting the development of the state of mechatronic process equipment, alerting 
personnel in case of emergency and other dangerous conditions, and, if necessary, updating or adjusting control 
programs. Provision is made for performing of some of the technical state maintenance functions by the mechatronic 
facility itself, 1.e., equipment self-service. The concept of “autonomous management of the technical condition” is 
formulated. The system structure and control functions are considered. It is noted that the implementation of the 
systems under consideration will significantly increase the efficiency of the equipment use. The performance of the 
autonomous control of the device or MPF in general is evaluated in accordance with ISO 13381-1: 2004. Based on this 
standard and the data presented earlier, a neural network structure is built to assess the autonomy of state management. 
The system training efficiency is considered taking into account the standard deviation of the network outputs from the 
target values of the training sample. 

Discussion and Conclusion. A list of the basic control functions at different levels of maintenance autonomy is 
presented: from alarm for failure prediction to complete maintenance autonomy without the direct involvement of an 


operator. 
Keywords: digital systems, autonomous maintenance, technical condition management, critical condition. 
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Introduction. For more than 50 years, serious developments have been carried out in the field of machine 
maintenance, which focus on methods and standards for maintaining an operational condition of the equipment. New 
devices and techniques are being developed for determining emergency conditions, identifying faults, and predicting a 
remaining life. The condition control and monitoring are essential for production mechatronic facilities of the industrial 
enterprises, such as machines with computer numerical control (CNC), forging and foundry machines, automated 
flexible transfer lines, etc. Equipment maintenance methods take various forms [1, 2], each of which has pros and cons 


that affect the level of control of the production condition of the equipment. Such control and providing the reliability of 
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machines are urgent scientific tasks. Numerous studies are devoted to this topic. An example is the works [3-14]. 
Russian and foreign scientists mainly consider the creation of automated monitoring and diagnostic systems. When 
using artificial intelligence methods for automation, processor modules are developed; they have acquired a general 
name of “e-maintenance systems”. 

Digital systems of electronic service are implemented by a collection of hardware-software components, which 
are designed to receive sensory information and assess the equipment status realtime. They provide predicting state 
changes in the production equipment and its components, generating messages on the need for certain actions on the 
part of experts, and, in some cases, changing the equipment operation parameters. The systems under consideration are 
used for personal diagnostics and monitoring of the technical condition of separate mechatronic objects. In addition, 
they are involved in the service system of a group of production machines. 

The focus of current research is the creation of programs and algorithms for service, and analysis of its results. 
Developments in the field of service quality and efficiency of the equipment utilization are aimed at this. In modern 
production sectors, such systems as Casip, Enigma, Icasame, RemoteDataSentinel, Intermor, INID, IPDSS, Proteus, 
MRPOS, WSDF, Telma, etc., are used for servicing [15—18]. 

It stands to mention multiplicity of issues on maintaining the reliability of the condition and arising faults, as 
well as the multifactority of the application of existing monitoring and diagnostic systems. In this regard, the tasks on 
the development of methods and techniques for the processing equipment maintenance, which provide decision-making 
under the conditions of data veracity and limitation, remain relevant. 

In recent years, the basic functionality of control, monitoring and diagnostics systems, tools, methods and 
technologies for servicing, forecasting and analyzing equipment status has changed. They are reoriented from the 
analysis of the technical condition to the prevention and detection of malfunctions at the early operational stages. Such 
control methods are being developed that allow minimizing the impact of emerging malfunctions on the performance 
and downtime of the equipment, on the number of failed parts and the frequency of failures. Ultimately, new 
approaches reduce economic losses. 

The urgent issue is the development of a new solution based on the intelligence of the system. It is about 
creating an automated autonomous service management and, as a consequence, about control over the autonomous 
technical condition of the equipment. 

Autonomous maintenance (AM) [19] is a comprehensive definition. In relation to machines, it implies such an 
automatic performance of service and maintenance functions, in which the equipment independently maintains its 
operability (self-service). If human intervention in the maintenance is wanted, the service functions remain hand-driven 
by an operator, that is, the problems are fixed promptly, without involving external repair services and until the end of 
the production run. 

The self-servicing of mechatronic process facilities (MPF) in the sense presented earlier is not yet an 
implemented solution, but a promising one. It is expected that the methods of artificial intelligence in the field of 
diagnostics and servicing will receive a level of knowledge equal to the level of staff. 

Materials and Methods. At the present stage of technology development, we consider the concept of MPF 
AM in the following interpretation. The purpose of the solution is maintainability based on an assessment of the state of 
the component and the MPF in general. For this, an independent automated control purposefully displays the status and 
required impacts on the MPF components [20]. 

As noted earlier, the basic approach to creating and developing autonomy in managing the operability of the 
processing equipment is the application of artificial intelligence and intelligent control methods under the conditions of 
uncertainty and fuzziness of the knowledge provided [21]. 

Research Results. Under the conditions of uncertainty, we consider the control procedure as the dependence 
of the real state on the goal state (successful performance): 


R*=AB, 
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where A and B are expert reports on the real and goal state of MPF. 

Presentation of the result of fuzzy inference Y depends on: 

— observed data_X obtained by the diagnostic system; 

— knowledge system 4A > B; 

— coefficient K, which is responsible for a possible measurement error. 

Therefore, the derivation of Y takes the form: 

Y=X x R*+K =X x (A >B)+ K, 

where Y = {Y,, Yo, ..., Ym} are parameters responsible for the equipment state management and for log output for the 
operator; X = {X,, Xo, ..., X,} are data on the input state of the equipment. 


Using this interpretation of the maintenance autonomy, we consider the structure and functions of the MPF 


AM (Fig. 1). 


pre eeerooroeeeeeeseee=e=~e; 





Device 1 Device 2 _ Device n 


Fig. 1. Structure of the control system of the equipment technical condition 





Here, signals from the diagnostic system sensors arrive at the control device not directly, but through the 
control unit of the technical state (TSCU) of the nodes. 

Nodes of the control unit receive data from the central node of the autonomous control of the process system 
(PSAC CN), which: 

— coordinates the operation and maintenance of all machine components; 

— processes data received from diagnostic systems; 

— on the basis of these data, analyzes the state, predicts the development of malfunctions and makes decisions 
on the required actions on smooth operation of the equipment. 

In the proposed approach to the generalized functions for the MPF technical condition management, control 
methods using signaling about the risk of failures, emergency situations, the need for a fuzzy device working band, etc. 
[22], are considered. 

An electronic service system (for example, e-Mind Machine) can provide information communication between 
the PSAC of an object through IS with equipment maintenance and repair units and enterprise management structures 
[23]. 

The autonomous control functions include decision-making on the system behavior in the event of a 
malfunction; consequently, the system should be trained during operation. Thus, diagnostic and monitoring methods are 
being improved, which provides eliminating previously unknown malfunctions. 

The analytical methods are complicated (or impossible) in adequate presentation of the state and changes of 
MPF. This is due to the following maintenance features of this equipment: 

— structural complexity of the mechatronic object consisting of a large number of nodes; 

— variety of factors affecting the reliability and performance, accuracy and productivity of machines; 

— wide range of properties of the materials from which the parts are made; 

— variety of requirements for functionality, parametric characteristics and performance of units, modules and 


devices in accordance with their purpose, accuracy characteristics, implemented production tasks. 
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The e-MindMachine concept provides management based on data on the real state of the facility and 
operational history. In this case, the development of faults for each unit is predicted. The application of the gamma- 
percentile life of reliable operation is proposed if statistical data are available or a border health band 1s used [22-24]. 

When developing PSAC of process facilities, it is possible to provide a ranking of autonomy levels depending 
on the maintenance unit used at the enterprise [25, 26]. It also depends on the criticality of the machine in the overall 
production. In this case, criticality determines significance of one or another of its units for the enterprise. The higher 
the criticality, the more economic, physical or temporary losses an enterprise suffers from the downtime of this unit. 
Such data is based on the equipment usage statistics, the number of failures, downtime, the cost of spare parts and 
manufactured parts, the number of work shifts, etc. 

When analyzing the criticality, it is also important to consider the future conditions of the machine, the degree 
of development of malfunctions, and the planning of the acquisition of spare parts. This is necessary to assess the 
residual life and to calculate the time before a possible failure under the conditions of uncertainty and stochasticity of 
data on the equipment status. Criticality assessment helps to find out the significance of the given equipment unit for the 
enterprise, and how much stop service and work resumption cost. 

To analyze maximum number of parameters that affect criticality of failures, we construct a system for the 
generalized assessment of the MPF technical condition criticality. For the output aggregate criterion of criticality, we 
use an ANFIS-based neuro-fuzzy assessment in MATLAB. As an example, we consider the failure criticality of the 


spindle assembly of a HAAS SMM vertical milling machine (Fig. 2). 


x22 0.5 











Fig. 2. Failure criticality assessment. Red vertical line is indicator 0.5: 
(a) equipment performance and accuracy; (b) costs of equipment repair and spare parts; (c) equipment downtime; (d) 
(b) number of equipment (unit) failures for a certain period; (e) time for fault identification and repair; 
(c) (f) Y=1 is full compliance to output parameters 

Fig. 2 shows that the input actions in the neuron-fuzzy system are responsible for all the key aspects which 
concern the equipment operation at the enterprise [27]. As a result, when logging-out, we obtain a generalized criterion 
of criticality Y whose value can vary from 0 to 1. Here, Y = 0, the input parameters do not meet the production 
requirements. According to the monitoring data, the parameters went beyond 0< Y< 1, but an expert can make a 
conclusion on the availability of equipment (for example, at Y= 0.75). Using this criticality assessment technique, when 
determining the operability of the same type of equipment, it is possible to indicate a response range for emerging 
malfunctions and create a new decision-making system under uncertainty. Thus, the most significant (critical) 
equipment receives maximum response from the operator and maintenance personnel for its operability. 

Fuzzy logic methods simplify the intelligent control system and evaluate immediately the state of an object. 
This is an autonomous management of the technical condition of the facility, so you need to understand how deep 
autonomy can be. A machine is considered together with an operator, which means that the equipment is automatically 


controlled by the actions of the operator as well. Management levels depend on the functions assigned to human. 
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Level 1. Control. It is determined by the presence or absence of an emergency malfunction, which causes 


equipment shutdown. 


Level 2. Diagnostics. The data analysis on the equipment condition is carried out by the built-in or mounted 


Sensors. 
Level 3. Monitoring. The current level of performance is determined, operational or tactical decisions are 
made. 
Level 4. The operator performs all basic functions without machine automation. They can be displayed as 
instructions. 


Level 5. The operator performs part of the functions that do not require deviation from the production tasks. 

Level 6. Autonomy of the MPF service without the intervention of an operator. 

If it is necessary to maintain the equipment operability at a high level, it is advisable to use monitoring (level 
3). The use of some automation at the 4th level seems to be a transition to the creation of intelligent service systems. 

The implementation of the systems under consideration increases significantly the equipment efficiency. The 
productivity of the autonomous control of a device or MPF is generally evaluated in accordance with ISO 13381-1: 


2004. Based on this standard and the data presented earlier, we construct the structure of a neural network for assessing 
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the autonomy of state control (Fig. 3). 





OUTPUT 






Fig. 3. Structure of neural network assessment of the state management autonomy’ 


Learning efficiency was evaluated with account of the standard error of the deviations of the network outputs 
from the target values of the training sample [20]. 

Discussion and Conclusion. The paper presents the basic principles and approaches that provide designing, 
creating and improving autonomous systems for management of the technical condition of the process facilities. The 
operation of such systems implies significant uncertainty and limited aprior and current data for decision-making. This 
necessitates the introduction of state-of-the-art technologies and methods of artificial intelligence. 

The authors’ view on the role, functions and construction of maintenance systems for MPF and other 
equipment is presented. The proposed approach is based on the application of digital information processing and 
decision-making technologies, which is due to the growing complexity of MPF automation. In each process system, the 
principles of openness, autonomy, and control of fine movements are implemented to one degree or another. This 
should be considered when designing workpart procedures and when implementing processes to maintain the 


equipment operability. 


' PPM — preventive and predictive maintenance. 
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Introduction. The artificial intelligence (AI) performance in data storage systems is considered. When working with 
data, the advantage of its use both in economic terms and for security is determined. The work objective is the introduc- 
tion of artificial intelligence in data storage systems. The key tasks involve the description of methods for data separa- 
tion, organization of its storage and counteraction to security threats. 

Materials and Methods. The data that should be fed into the drives is divided into parts so that it can be restored without 
one of the parts. This is necessary to be able to access and recover information in the event of a software or hardware 
failure. 

Results. The AI performance under detecting security threats is considered. Since the model implies the interaction of 
users with data, it was found out how the data access control is carried out and the keys are stored. 

Discussion and Conclusions. The use of AI in organizing a data warehouse will speed up the system. Artificial intelli- 
gence with built-in machine-learning algorithms will provide responding to a situation that affects the state of the sys- 
tem. Analysis of the state of the drives will avoid a possible hardware or software failure. Minimization of the human 
factor in the system operation contributes to the improvement of its work. 


Keywords: artificial intelligence, threshold separation, threat, machine learning, storage organization, dynamic change, 
key, encryption, backup, attack, hash. 


For citation: V. V. Zhilin, O. A. Safar’ yan. Artificial intelligence in data storage systems. Vestnik of DSTU, 2020, vol. 20, 
no. 2, pp. 196—200. https://doi.org/10.23947/1992-5980-2020-20-2-196-200 


Introduction. At the present stage of development of information and telecommunication technologies, 





a per- 
son comes across various kinds of information. The data storage in a safe place is an urgent issue. The science that deals 
with the storage of information in an encrypted form is cryptography. It studies how to hide data and provide its confi- 
dentiality’. 

The information is static in all data storage implementations. That is, the data will be where it was for the first 
time. Information can get to another disk or another sector only after removing it from its original location. However, 
this fact may be regarded as weakness of existing algorithms. 

To increase the protection effectiveness, it 1s recommended to share data [1] in order to store their parts. This 
provides organizing fault-tolerant storages for information recovery algorithms in the event of any fault, failure of one 
of the drives, as well as with loss of one or more parts of the data. 

If an attacker has access to the disk, he will be able to obtain the necessary information. In case of static stor- 
age, if an illegitimate user gains access to several such disks, he can recreate information. Separated data storage in- 
creases information security. In addition, in large storages, as a rule, hard disk drives are used, whose running speed 
(reading and writing operations) is significantly lower than that of solid state disks (SSD) and flash drives’. In this re- 
gard, you have to choose between the volume used and the operation speed. 

The work objective is to describe the performance of artificial intelligence (AI) in data storage systems. The 
task is to describe the AI operation algorithms when performing the following operations: 


¢ data recording on drives; 


'Ryabko BYa, Fionov AN. Kriptografiya v informatsionnom mire [Cryptography in the Information World]: Goryachaya liniya-Telekom; 2018. 305 
p URL: https://www.techbook.ru/book.php?id_book=1001 (accessed 04.05.2020). 

*Sravnenie SSD i HDD diskov v real'nykh usloviyakh ispol'zovaniya [Comparison of SSD and HDD drives under real-time conditions of usage]. 
Habr. URL: https://habr.com/ru/post /394135/ (accessed 04.05.2020). 
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¢ referral from users; 

¢ data warehouse analysis; 

¢ in case of information security threats. 

Artificial intelligence impact on the storage system quality. In modern devices, such as smartphones and 
computers, developers give special consideration to the implementation of artificial intelligence. In these technologies, 
machine learning algorithms are implemented, which increases the operation speed and reduces the response time to 
conducting frequently repeated information. It should be considered how the introduction of artificial intelligence and 
machine learning into storage systems can significantly improve the quality of their operation [2]. 

Currently, large companies and organizations are starting to implement AI in their data warehouses. A wide- 
spread adoption of machine learning and artificial intelligence technologies helps to improve the operating quality at the 
management level. This facilitates the work of network and data warehouse administrators through constantly diagnos- 
ing the causes of traffic overload and reduction. This will allow them to identify potentially vulnerable segments of the 
model beforehand. 

Artificial intelligence involves the use of integrated deep learning algorithms that can predict the state of the 
entire system and respond quickly to possible changes. This will significantly reduce the cost of eliminating the conse- 
quences caused by equipment failure. In addition, the introduction of artificial intelligence in the organization of fault- 
tolerant storages will allow for their automation [3]. This implies an analysis of the system state and processing of in- 
coming data in dynamic mode. 

Artificial intelligence and machine learning provide minimizing the likelihood of data loss. Together with re- 
dundant arrays of independent disks, such a system increases the availability and speed of overcoming forced downtime 
thanks to intelligent data recovery, a backup strategy and the transfer of necessary data [4]. 

Data distribution by input parameters. To store data in a shared form, threshold separation methods can be 
applied. In classical algorithms, the in parameters are static values, which is a major disadvantage. Gaining access to 
one data through selecting input parameters threatens safety of all other data in the system. 

Thus, the greatest preference in terms of information security is given to the algorithm that uses various input 
parameters for threshold sharing. These parameters can be randomly generated according to a certain algorithm or de- 
pend directly on the input data, which will be analyzed accordingly to select the most preferred parameters. The gen- 
eration of such parameters will be handled directly by artificial intelligence. At the same time, it should consider the 
number of users to whom this information is available. Threshold sharing parameters and the location of the first share 
are stored in the database. This database is encrypted on the keys of specific users, which corresponds to the application 
of an asymmetric cryptoalgorithm. 

Parameters of the proposed algorithm for data reliability. When distributing data among drives, artificial 
intelligence uses the following parameters: drive speed, availability, free volume and reliability index. At the same time, 
Al uses the data value parameters, size and frequency of accesses to it for calculation. 

The algorithm for the distribution of data shares across drives is based on the calculation of drive coefficients. 
To determine the speed of the drive S, it is required to find the arithmetic average of the writing speed s,,,, and read rate 


San 
S3antSyr 
2 
Drive accessibility A has 3 levels: 0 — unavailable, 1 — frequent situation of drive unavailability, 2 - rare sit- 


S= 


uation of drive unavailability, 3 — anytime available. The exponential coefficient of the drive can be expressed by the 
formula: 
K, =S:-A-V-R, 

where V is disk capacity; R is reliability, it is ranked from | to 10 points; A is the number of file accesses in a given peri- 
od of time. 

The significance factor of the file Ky depends on the value level (1 — low, 2 — medium, 3 — high). Based on 
the known attributes, this parameter can be calculated from the formula: 

Ky =S°V-A. 

Consider the technique of selecting a storage drive for the remaining files. To determine the priority of the 

drive, we use the file number n, in the sorted table in descending order of the parameter Ky. Then the disk number 


N, participating in the sample for storing the first share of the data is found from the formula: 
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d . Ny 
N, = round (—) mod no, 


where d is the number of disks used for storage; fis the number of files in the system; nz is the number of drives; round 
is rounding to the nearest integer. 

In the same way, drive NV is determined, which may be selected as a drive for storing the first share: 

Kyirst = Kg Ka max, 
where Kg is the split file ratio; Kg may 1s maximum drive coefficient. 

Next, the drive, whose difference in the coefficients Ky;,,, and K,, is the smallest, is selected as N>. To store 
the first share, one of the found disks NV, and N> is randomly selected. 

Advantages of the considered algorithm of data distribution among storage devices. Classically, disks 
combined into one data array (RA/D-array) operate according to a static algorithm. That is, when the first share of data 
is detected, a search procedure for the second share can be performed. Each time, upon detection of the next share, the 
probability of determining a specific distribution algorithm increases. In the proposed algorithm, the distribution 1s 
based on information data and drives. In other words, the task of finding all the data shares is NP-complete, that is, it 
cannot be solved in polynomial time. It follows that the use of this algorithm increases the reliability of the information 
Storage system. 

However, users do not create a key; it is generated automatically and stored on the devices of the users them- 
selves. Thus, when trying to decrypt data, a specific key must be used, otherwise this operation fails. The location of the 
keys changes dynamically over time. Such changes can occur at certain hours, or at specific time intervals [5]. Artificial 
intelligence is responsible for the key location with account of the availability of all devices connected to the system. 
This complicates the work of an attacker whose goal is to access information stored on drives [6]. 

Since the information in the repository is operated by users, it should be determined who of them is allowed to 
conduct operations with data, in other words, how to differentiate the data access. 

Access comparison is performed according to the hash-sum table’. If a user has access to data, then they are re- 
stored. If there is no entry in the database on granting access rights, the entity requesting access to the object is denied. 

In addition, an attempt to gain access will be recorded in the log file, and the subject who owns the information 
stored in the system will be informed about the access attempt by an illegal user or subject. In this case, artificial intelli- 
gence analyzes the actions of each user to make decision in the event of a specific situation. For example, it determines 
whether the query is invalid, or it has any purpose [2]. 

Fig. | shows an example of granting access to data recovery and access denial with recording the event in a log 


file and the user notice. 





at, 
es 
Fig. 1. Granting and denial of access to data recovery 


Consider the access control function. Since the primary subjects in the system are users, the issue of providing 
access to data is particularly topical. Full control over user operations is very important. The decision to provide access 
to data is made by a distribution device with support of artificial intelligence. In this case, the use of AI will increase the 


' Legkii sposob poddelki kontrol'noi summy s pomoshch'yu kollizii [Easy way to fake checksum using collisions]. xakep.ru. Available from: 
https://xakep.ru /2012/11/22/light-fake-checksum/ (accessed 04.05.2020). 
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user’s data handling speed. This is achieved through transforming the most commonly used information, which avoids 
waiting for the response time to the sent data recovery request. 

We consider dynamic change of storage. In certain cases, artificial intelligence controls the transfer of data and 
is responsible for data backup and restore including the time-based ones. In this case, the data migrates from one disk to 
another, and the information in the database also changes. This minimizes the likelihood of predicting the system state 
at a specific time. Since all keys are generated on the ground of previous sequences using hash sums, the task of finding 
a key designed to decrypt a database is an NP-complete problem [7]. In this case, artificial intelligence analyzes the 
state of drives, as well as the system as a whole. When a malfunction is detected, data is transferred from those drives 
on which failures are noticed. This minimizes the likelihood of information loss. 

When a threat is detected, artificial intelligence analyzes the attack to determine its target. If the final goal is 
the storage object, then the data from it is transferred to another drive. Thus, in the event of a successful attack, AI will 
take effective actions to hide information [8]. 

As users work with data, the human factor plays an important role. Therefore, the system is designed so that 
only users with access can work with the data. There is an analogy with the mandatory access control model. An excep- 
tion is that the distribution of access is controlled by the user who has created the information in the system and is its 
owner. 

Conclusion. Thus, the use of artificial intelligence increases the system speed under organizing a data ware- 
house. Restricting user access to the algorithm of the system operation improves information security. Artificial intelli- 
gence with built-in machine learning algorithms provides a quick response to any situation affecting the state of the 
system. Analysis of the state of the drives avoids a possible hardware or software failure. Minimization of the human 
factor in the system operation helps to improve its functioning and a deeper analysis of user queries. In addition, the 
collection of information on possible attacks enables maintaining the system security at a proper level. 
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